Équation différentielle de Lamé et fonctions de Lamé 


L'équation de Lamé est un cas particulier de l'équation de Heun, elle s'écrit : 


Sous sa forme algébrique : 
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Sous sa forme Jacobienne : 


. D Ze 2 . . * 
en considérant z = sn“(u,k), il vient : y"'(u)+(h=v(v + Desn (u,k))y(u)= 0 


Sous la forme trigonométrique faisant intervenir l'angle o de l'amplitude Jacobienne : 


Q= -~ am(u,k) = Соѕ(ф) = sn(u,k) 
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Sous la forme donnée par C.G . Lambe (1952 < Lamé-Wangerin functions > ) : 


Dans l'article de 1952 de C.G. Lambe, « Lamé-Wangerin functions », les fonctions sont 
improprement dénommées « Lamé-Wangerin >x, alors qu'il s'agit des fonctions actuellement 
dénommées < fonctions algébriques de Lamé > ayant pour propriété d'être de période 8K et pour 
certaines valeurs du paramètre h de développement fini (la raison pour laquelle on les appelle 
algébriques) : 
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Avec beaucoup d'algébre et gráce à Mathematica, l'équation de Lamé en x et u(x) s'écrit : 


x(x 1x а) Ë ier 2(a H 1)x H 22 H DEI Dx À v?(a H 1)u(x)=0 


dx? 


En posant v=u+1/2, il vient : 
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Cette dernière équation sert pour étudier les solutions polynomiales en z lorsque n est un entier, et 
dans ce cas les solutions de l'équation de Lamé sont dítes algébriques. 


Solutions périodiques 2K et 4K de l'équation de Lamé 


De la forme Jacobienne de l'équation de Lamé, il apparaît clairement qu'il s'agit d'une équation 
différentielle du second degré dont les termes sont périodiques. Il s'agit dont d'une équation de 
Hill, mais au lieu d'avoir des termes sinusoïdales de période 2л ou x, il s'agit d'un terme sn°(u,k) 
impliquant une fonction de Jacobi de période 4K, soit de période 2K puisque le terme est au carré. 
Aussi paraít-il naturel dans une premiére approche de rechercher des fonctions solutions de méme 
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Ces fonctions périodiques sont d'une importance considérable car elles sont souvent parfaitement 
adaptées à la géométrie induite par les systèmes de coordonnées orthogonales utilisés qui 
présentent la périodicité 4K ou 2K dans certaines de leur variable. 


On peut à profit établir une analogie avec soit les solutions périodiques de l'équation différentielles 
des vibrations, soit avec l'équation de Mathieu. Cette analogie permet de clarifier (sans pour 
autant expliquer et justifier) la construction des solutions périodiques de l'équation de Lamé. 


Il faut vraiment se reporter aux articles fondateurs dans la construction des solutions périodiques 
pour comprendre le « comment du pourquoi », étape je pense indispensable pour toute utilisation 
adéquate des fonctions de Lamé dans la résolution d'un probléme aux limites, soit : 1940, E.L.Ince, 
«М The Periodic Lamé Functions >x, 1940, E.L.Ince, < VIl.Further Investigations into the Periodic 
Lamé Functions », < A.Erdélyi.H.Bateman, HIGHER TRANSCENDENTAL FUNCTIONS VOL III, Lamé- 
Heun Functions, chapitre XV, section 15.5.1 ». 


Il est bien évident que l'équation de Lamé possède des solutions d'un tout autre type, soit non 
périodiques, soit de période un multiple entier quelconque de 2K. Il est tout aussi évident qu'en tant 
qu'équation différentielle du second ordre, elle admet une seconde solution linéairement 
indépendante. Dans ce qui suit ces secondes solutions ne sont pas décrites. Un des moyens les plus 
simples de construire les secondes solutions est d'utiliser le développement en série de fonctions de 
Legendre de premiére espéce, puis de rechercher un développement similaire à l'aide des fonctions 
de Legendre de deuxiéme espéce pour les obtenir (voir 1948, Erdelyi « XXVII. Expansions of Lamé 
Functions into Series of Legendre Functions » et 1966 B.D.Sleeman, « The expansion of Lamé 
functions into series of associated Legendre functions of the second kind »). 


Dans ce qui suit nous établissons trois tableaux des analogies de construction des solutions de 
l'équation de Lamé, de l'équation des ondes et de l'équation de Mathieu, sur une méme ligne (en 
lecture) : 


Équation des ondes 
y" (z)*hy(z)20 : un paramètre h 
Notation Valeur Parité autour de | Parité autour de | Conditions aux | Période Translation л | Translation r2 

caractéristique h | z=0 2=п/2 limites sur 

[0,172] 

Cos(2nz) h=(2n)° y(-z)=y(z) y(rez)-v(z) у'(0)=у'(л/2)=0 |n у(2+л)=у(2) у(2+п/2)=у(2-п/2) 
51п((2п+1)2) h=(2n+1)° y(-z)=-y(z) у(т-2)=у(2) y(0)-y'(1/2-0 |27 y(zer)--yz) ` | у(2+п/2)=-у(2-л/2) 
Cos((2n+1)z) h=(2n+1)? y(-z)=y(z) у(л-2)=-у(2) у'(0)=у(л/2)=0 |27 у(2+п)=-у(2) | у(2+п/2)=-у(2-л/2) 
Sin((2n*2)z) h=(2n+2)? y(-z)=-y(z) у(л-2)=-у(2) у(0)=у(л/2)=0 |n у(2+л)=у(2) у(2+п/2)=у(2-п/2) 


Équation de Mathieu 


y" (z)*(h-2qCos(2z))y(z)-0 : deux paramètres h,q 


Notation Valeur Parité autour de | Parité autour de | Conditions aux | Période Translation rt | Translation л/2 
caractéristique | z=0 2=п/2 limites sur 
h=f(q) [0,7/2] 
се›.(9,2) hza»(q) y(-z)=y(z) y(rt-z)=y(z) y(0)-y(1/2)-0 n у(2+л)=у(2) |у(2+п/2)=у(2-п/2) 
Sexa[q,z) h=bxa(q) y(-z)=-y(z) y(rez)-v(z) y(0)-y(1/2-0 |2m у(2+л)=-у(2) | у(2+л7/2)=-у(2-п/2) 
Cexa(q,z) h=axa(q) y(-z)=y(z) y(rt-z)=-y(z) yv'(0O)=y(1/2)=0__|2n у(2+л)=-у(2) | y(z+1t/2)=-y(z-11/2) 
sex.o[q,z) h=b22(q) y(-z)=-y(z) у(т-2)=-у(2) у(0)=у(л/2)=0 |n у(2+л)=у(2) |у(2+п/2)=у(2-п/2) 


Équation de Lamé 
y''(z)+(h-v(v+1)k2sn2(z,k))!y(z)=0: trois paramètres h,v,k?, changement n/2->K et intervertion de Ес,2"*(2) et Es,2"**(z) dans la notation NIST 
Notation Valeur Parité Parité autour de | Conditions aux | Conditions aux | Période | Translation 2K | Translation К 
Nist/Ince caractéristiqu | autour de | zzK-K(K?) limites sur | limites sur 

e h=f(v,k?) z=0 [0,K(k2)] [0,2К(К?)] 
Ec/^(z/Ec/^(z) | hza/"(K) y(-z)=y(z) | y(2K-z)=y(z) y'(0)-y'(K)-0 y'(0)=y'(2K)=0 | 2K y(z+2K)=y(z) | y(z+K)=y(z-K) 
Ec2n"1(z)/Es2""1(z) | һ=а, (K2) y(-z)=-y(z) | y(2K-z)=y(z) y(0)=y'(K)=0 y(0)=y(2K)=0 | 4K y(z+2K)=-y(z) | y(z+K)=-y(z-K) 
Ез" (2)/Ес D'al | h=b,?" (k?) y(-z)=y(z) | y(2K-z)=-y(z) у'(0)=у(К)=0 y'(0)=y'(2K)=0 | 4K y(z+2K)=-y(z) | y(z+K)=-y(z-K) 
Es/""(z/Es,""(z) | һ=Ь, (К?) y(-z)=-y(z) | у(2К-2)=-у(2) y(0)=y(K)=0 y(0)=y(2K)=0  |2K y(z+2K)=y(z) | y(z+K)=y(z-K) 


Orthogonalité et Normalisation des fonctions périodiques 2K et 4K 


Il est assez facile d'établir les formules d'orthogonalité des fonctions périodiques de Lamé pour une 
même valeur caractéristique h : 


асва (x, k? )Ec2" (x, k: )= 0 pour m= p fas gem (x, k? )Ec2#" (x, e) 0 pour т=р 
0 0 
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La normalisation dépend des articles, mais je prendrais la normalisation la plus communément 
utilisée, celle donnée dans < NIST Handbook of Mathematical Functions, chapitre 29, Lamé- 
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La normalisation proposée par Ince dans « 1940, E.L.Ince, VII.Further Investigations into the 
Periodic Lamé Functions » semble être la suivante : 
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Le complément de normalisation fixe des valeurs positives aux fonctions ou leur dérivée première : 
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Algorithme de construction des solutions périodiques 2K et 4K de l'équation de Lamé par le 
développement trigonométrique 


Tout comme pour les fonctions de Lamé, la construction des fonctions périodiques 2K et 4K peut se 


présenter sous la forme de la recherche des valeurs propres et vecteurs propres d'une matrice tri- 
diagonales. Le développement en série le plus commode est celui en série de Fourier de la variable : 


| = > am(z,k) <> Cos(@)= sn(z,k)= ар = —dn(z, k) 


Le schéma de construction des valeurs propres et vecteurs propres d'une matrice tri-diagonale est 
toujours le même, et se présente génériquement comme suit : 


Pic. uS. ә 2а, „б ышы | 
п Ba 0 » 0r. us 
Quum. cue. 44 0 = 0.3: 
(M, - HI,)X, =0 à à 
avc Hz2h-v(v41) M, =|" 2 буз s à " Ress) 
Sc er Улу "un Dua y E 
h- LC 0 0 ess 0 у, bya me 
0 0 Yx B, 
Ec?"(z,k2)= S x. „(h „(h » JCos(2po) ou dn(z,k)Y х, (8, JCos(2po) 
p=0 =й 
| pz H, +v(v +1) | pco : pco 
H, valeurs propres — h, 2 Ес") Xy, „(A, ‚)Соз((2р+1)р) ou dal HI. Tos (h, JCos((2 p +1)o) 
X, (h,) vecteurs propres de M, > a d d 
melt? N...) Ез", k dÉ S x. „(h „(h,)Sin(2p+1)p) ou dn(z,k) DG „(h JN. )Sin((2p+1)p) 
p=0 p=0 
p=% 
Es?" (2, k?)= Xs; „(h E Jsin((2 p + 2}p) ou dnl (z, k) js, MUS }Sin((2 p + 2)p) 
p=0 p=0 
Yo=0 x, #0 quelconque fixé par la normalisation 
, Box, + G x, = H x, 
Récurrence 
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Les récurrences sont celles décrites dans l'article fondateur de 1940 de E.L.Ince < VII—Further 
investigations into the periodic lamé functions>. Dans cet article la variable d'amplitude prise est la 
suivante : v = am(z, k) = Sin(v)= sn(z,k) et Cos(v)= en(z, k) et dv = dn(z,k): Puis ce sont les 
récurrences avec la variable d'amplitude ф dans l'ouvrage < A.Erdélyi.H.Bateman, HIGHER 
TRANSCENDENTAL FUNCTIONS VOL Ill, Lamé-Heun Functions, chapitre ХМ, section 15.5.1», ainsi 
que celles de « NIST Handbook of Mathematical Functions, chapitre 29, Lamé-Functions ». 


Voici un tablequ exprimant les huit types de développements trigonométriques possibles des 
fonctions de Lamé périodiques 2K et 4K (deux par fonctions), avec les cas pour lesquels le 
développement est fini, soit qui donne lieu aux polynômes de Lamé correspondants : 


Fonctions de | Type de développement et formule dans le | Récurrence Normalisation Cas d'un développement fini : 
Lamé livre « A.Erdélyi.H.Bateman, HIGHER polynóme de Lamé et type et 
périodiques TRANSCENDENTAL FUNCTIONS VOL Ill, notation d'Arscott, type de NIST 
Notation NIST Lamé-Heun Functions, chapitre XV, section 
15.5.1», Repris dans « NIST Handbook of 
Mathematical Functions, chapitre 29, 
Lamé-Functions ». 
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B, SC sdE? (2,02) me... n-1 
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Formule (21) 
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Formule (24) 
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"à se dÉ vk KE (2p+2)0) 
Formule (22) px] 
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Formule (23) AN 


II est facile de se convaincre de l'équivalence entre la notation d'lnce (j'y ai rajouter d lorsque le 
développement alternatif commençant par dn(z,k) est proposé) et celle de NIST, à travers les 8 
développements présentés comme suit : 


Fonctions de | Type de développement dans le livre | Fonctions de Lamé | Type de développement tel que dans l'article d'E.L.Ince 
Lamé « A.Erdélyi.H.Bateman, HIGHER | périodiques Notation Ince, « Vil—Further investigations into the periodic lamé 
périodiques TRANSCENDENTAL FUNCTIONS VOL functions» ; „ = am(z,k) 
Notation NIST | 11, Lamé-Heun Functions, chapitre 
XV, section 15.5.1», Repris dans 
« NIST Handbook of Mathematical 
Functions, chapitre | 29, | Lamé- 
Functions >. ç = 5 _ amlz,k)= 254 
2 2 pem 2 2 
Sch) ER EE] LEN NE 
p=! p= 
Е) вов) Gies capo] ckt Е. jj 
р=\ р=\ 
2m+l 2 p=% 2m4l 2 E 
SKI | кые) S 4,„,cos(2p+1jp) E hël ве" Ч(ые)- S 4, CU sip) 
p=0 p=0 
E oil k) dE au | k) 2 2 V 
c, Ú FLH dr (ЫЎ Cp, со{(0р+1р) 5, Ú, Ex" (P )- di) C, „(У (рч) 
ri р=0 
Es?" (zi?) Es?" (2.k2)= NEL 2р+1)р) E" Gp) Gel )= ВС!) )'Cos((2p +1)u) 
2m+1 2 pa 2m+l 2 pem 
Es," (2,4 ) ES" (e k) dale KÈ D. Sisi dEc;" LE) Ee (z )- [z k) D... C Cos((2p+1)u) 
p=0 p=0 
2m+2 2 2m+2 2 p=% 
Es" (2.2) в) 5, sepa iter] DEENEN 
p=0 p=0 
2m+2 2 ps 2m+2(_ 12 
Es," (ск E (Je dalz, NY D, .Sin(2p+2h) dEs, [t | E (e )= dn(z,k NE "a 1" Sin((2p+2)u) 
p р) 


Remarque:il est clair que tous les coefficients des développements пе diffèrent que d'un signe près 
(-1}. D'autre part dans la notation d'Ince, les propriétés de parité suivant z sont plus cohérente 
avec l'emploi de la fonction sinusoidale de parité usuelle. De ce point de vue la notation d'Ince est 
plus cohérente que celle utilisée par NIST (interchange des Ec et Es d'ordre impaire) et rejoint celle 
des fonctions de Mathieu (Ec, Es pour ce, se). 


Algorithme de construction des solutions périodiques 2K et 4K de l'équation de Lamé par le 
développement en fonctions elliptiques de Jacobi 


Les 8 fonctions de Lamé périodiques 2K et 4K sont essentiellement les mêmes que les précédentes. 
Il s'agit simplement d'un développement et d'un facteur de normalisation différent. On peut par 
exemple choisir de poser le premier coefficient du développement dans la sommation comme étant 
égal à 1. Cela entraïne bien évidemment une fonctions de Lamé qui diffère uniquement d'un 
facteur de normalisation. La récurrence est celle décrite dans l'article fondateur de 1940 de E.L.Ince 
« V.—The periodic lamé functions >. 


Le schéma de construction se présente génériquement comme la recherche des valeurs propres et 
vecteurs propres d'une matrice tri-diagonale : 


h, valeurs propres 


m 


x,h,) vecteurs propresde M, > 


m € (012... N...) 


[ & 0 Ss zm (( E 
п D o 0 X (ue 2 
0. uu ux ы; 0 L 0 … 


0 ууз Bra б; O moa 
0 у В Oya 0 … 


r=o% 


r=0 


г=%0 


г=0 


r=0 


r=0 


H 
8 


Yo=0 x,=0 


x =l < normalisation 


Récurrence 4 Byxy t oux, =h x, 


У.Х, + Bx, +a,x 


= h x, е „Ху + В,х, + VX, e h x, 


r rad 


Х, = Е 


r=% 


Ee” (zk )= X, n, )sn* (z,k) ou dEc?" [z k) )= ап(г, ЮУ x, „(h„)sn* (z,k) 


r=0 


r=% 


Evert X, „(h„)sn*"(z,k) ou dEs?""(z,k*)= dn(z,k) x 


r 


"" 


Ec?" (zi? )- cn(z,k) x, Un, ui (z,k) ou dEc?""(z,k2)=cn(z,k)dn(z, k) 


r,m 


r=% 


r=0 


> X (h, Jn” (z, k) 


(h, ui? (z, k) 


Fonctions de | Fonctions Type de développement | Récurrence Récurrence Cas d'un développement | Récurrence du 
Lamé de Lamé | / Normalisation forme 1 forme 2 fini : polynôme de Lamé | polynôme de Lamé 
périodiques | périodiques et type, notation et selon Arscott: ^ 
Notation Notation développement change de signe et 
Ince NIST/Erdelyi d'Arscott assignation de y. ^, 
ne change pas de 
signe 
к!) Fckr Ec?" [z, p dÉ + sn” (2, k) 24, +h4, =0 a, =-(2r+1}2r +2) v=2n>0 nh a, =(2r+1)2r+2) 
v = (r2), B.=4r (+) ре Sco has B, Arte) 
Normalisation д =1 (4206), – 22-2442 (uwen: j 7,7 uar 22r - 1€ 
-(v-2r+2)v+2r-1)k*4, , 20 
gees uE e) me un) 
ие) S CI Asn” (esk) 
r=0 
2m 2 2m 2 5 Е Le š 26,4 (i-e) =0 а, =-(2r+1)2r+2) ES 2 1 m a, =(2r+1)2r+2) 
deit Ec"{z,k авс) У C (д) | 26 ° Е и V=Zn+l Fa "H { 
v e (2r+1Y2r+2)C,,, + B, 4r +(2r +1} k? В 4 B, =4r"+(2r Je 
Normalisation С =1 а-а Dram y, cerea pones (Arscott), Type 4, гузора днее 
-(v-2r «fv EE 
(v-2 +v Gs d ЫР) тєл 
a£; EE} de Cr Cm) 
r=) 
mif, 12 mf, 12 "m m , 24, +(h-1)4, -0 а, =-(2r +1X2r +2) = = а, =(2r +1)2r +2) 
Ec" ke) Es, LE) Ec" ie )- o YA si" (z,k) MA u -—— Wc Зу, z 22 
` GC" (r+1f2r+2)4,,,+ B, =(2r +1} +4r°k Type 3 (Arscott), Type 3 | |? =(2r +1) +4r°k 
Nomnslisetion E И) E y, =(v-2r +1)v +2r)k? (NIST) : y, =-(2n-2r+2)Y2n+2r +1) 
-(v-2r+1)v+2r)t24, =0 c В) me Lal 
dir m ne, Dy -I Aal 
1) 
dE Ec" (z 73 Eg") aic" C Jen Cm 26, «(i -1- jo, =0 a Een ca у=2л>0 ra=- о) 
| 2 (2r +1)(2r+2)C2,.2 + B, - Gray (+) Type 5 (Arscott), Type 7 | |^ =(2r+)) (+z?) 
Normalisation c =1 (а-а (ее), = | xov 270 | (NIST) : y, 7-Qn-2ran + 2r + 
BET EE aE") wel, A 
n rtl 
«dE 28у HJ c sa? (2,4) 
= 
Eget e) Ec" (zi?) к sig) sn sn[z EB sn? А3) 2x38, з e gett EE HR De 
EH HH в, =) Туре 2 (Arscott) Туре 2 в, cer afl) 
Normalisation B =1 р-н) (se WW y, 2 (v-2r «Xv +2r)k (NIST) : y, =(2n-2r «22n 2r +1)k 
-(v-2r «Tv e 2r) B, , -0 E p тє (б.п) 
st emet] НЗ, 0) 
rz 
dE HIE Ü E 2т+1 k? E CREE 2x3D. (0 j= di =-(2r+22r +3) v=n>0 r zn-l a, =(2r +2)(2r +3) 
A Ú, C, (2, dE" ишы], NU Tz mx T ds 
! "d Dean , БАЎ +(2r+2J k? B, = Qr +(2r+2Y k 
Me 29 a | Type 7 (Arscott), Type 6 „2 
Normalisation D, =1 +h-( 2r+1) -(2r+2fk vk, — Se 2rv +27 +1) (NIST) : y, = -Qn-2r(2n «2r +1 
-v ponen D, -0 sdE! (zk) ae, al 
sii рар Јр, Cl) 
m 
2m+2(_ 12 2m42 ERT ER 4-6), = r = AMETE 
Eck El Es," [zk CAR а), d) Zeg P " ба E 5 deb o d SS Re 
2 (+22 +3)В„.,+ В.= 0+2) sre | Туре 6 (Arscott) Type 5 | 18 =(2r+2} «Qr e 
Normalisation B =1 «(i-r -orfe a- т, =(v-2r)v+2r+1)k? (NIST) : у, =-(2n-2rY2n+2r +1} 
-(v-2rYv+2r+1k2B, =0 at iz (zk 1 mel... n) 
rtl 
s ET Eo, Ui It, st) 
dE" (zk?) Eg) dg" P -ad ul so, si ) (2204-40, =0 a, =(2r+2J2r +3) v=2n+1>1 Ge =n-1 | (a, =(2r+2Y2r+3) 
(2r+2)2r+3)D,,.3 + DË Type 8 (Arscott) Type 8 (B - r2) +) 
Normalisation; y ed REECH BEE EE , у, =-(n-2rY2n+2r+3)e 
-(v-2r-1fv+2r+2J*D, =0 alt. Lu ach, n) 
soll k n dilz Eje li) ife D, sit) 


Exemple des premiëres fonctions de Lamé de périodes 2K et 4K développées en fonctions 
elliptiques de Jacobi 


Les premières fonctions de Lamé telles que développées par Ince sont aisément calculables car les 
valeurs propres h ont une expression relativement simple. Il s'agit évidemment également des 
expressions des polynómes de Lamé puisque le développement est fini. 


v=0 m-0 h-0 Ect(z,k2)=1 
v=l m=0 h=k dEc(z,k*)=dn(z,k) 
v-2 m-0 Acht Vies) кейш? ]=1-(гый ie keo 
v-2 m-2 һе) Ec?(z,k2)=1-(1+k2+V1-k? +K" Je pi 
v=3 m-0 й=2+5К°-21-К°+4&*% dEc(z,k°)=dn(2,k ha - ie ea bn?(z,k) 
v=3 m-2 h=2+5k «241- E? +48 = dEc?(z,k*)=dn(z,k) Aa, AC + 4k" bn?(z,k 
у=1 т=1 Л=1 Ecl[z,k2)= en(z,k) 
v=2 т=1 h=l+É dEc\(2,k°)= dn(z, k)en(z, k) 
v=3 т=1 h=5+2k*-2V4-k*+k* Ecl[z,k2)= cn(z,k (к) br JA- + bert AN 
v=3 m=3 h=5+2k «244- Kr В n(z, Qi -(244244/4-1 E wo) 
v=4 m=l h-545k 244 +48* _ dEc\(z,k?)= dn(z, k)en(z, D DEI VA + 4k" |a (zt), 
v=4 m-3 h-545k «244? +4k* dEc)[z 2 )= d nf, Hjen(z, E) - 278 a+ + 4k" Jsn?(2,k 


Е 


v=l m=1 h=1+k Es![z,k2)= sn(z,k) 
v=2 m=l h=l+4É dEs!\(z,k*)=dn(z,k)sn(z, k) 


v=3 m=l h=5+5k 244 T? + 4k взел) me) Lee ae ET a t) 


v=3 m=3 h=5+5k «24A- T) +4К* взе) mie ki Ab är ea m c Gt) 


v=4 т=1 h-5410£ -2J4-9E? +94 авзо) die tm) 1-2(2+ 382 = aco em ra) 


v=4 m=3 h=5+104°+2V4-94° +94 ate?) dr nl t) i +32 + V-94" e nea) 


v22 m=2 h=4+k Es?(z,k?)= cn(z,k)sn(z,k) 
v=3 m=2 h=4(1+k*) dEs2[;,k2)= an(z,k)en(z,k)sn(z,k) 


v=4 m-2 h=10+5k3*-2V9-9k*+4k* Es2[; k2)= EE EE V9 — 
v=4 m=4 h=10+5k3+2V9-9k3+4k*  Est(z,k?)= ам) -163+242+49- oC AC kt z,k 

v=5 m-2 h-2l5«5E 3M - i£ +4") dEs?(z,k?)=dn(z,k)en(z,k)sn(z,k)1 Ja А +k ka? (z,k) 
v=5 m-4 h-2l5«5E «341-8 «k*). dEst(z,k?)=dn(z,k)en(z,k)sn(z,k)1-(1+k2+vV1-k2+k* bn? (2,4) 


Au delà les valeurs de h sont des racines de polynômes de degré 3 et plus, qui n'ont pas 
directement d'expressions algébriques. 


Constructions des solutions de l'équation de Lamé de période purement imaginaire 2iK'et 4iK' 
Il s'agit ici de construire les solutions de période imaginaire de l'équation de Lamé : 
y"(z)+ (h -v(v + D) sn? (z, DÉI — 
telles que b (z 2iK())- y(z) k'=\1-k2. Pour cela la construction est assez simple, il suffit 


ou y(z + AiK (k')) = y(z) 
d'effectuer un changement de variable dans l'argument z sachant que : 


ТАЗЕ О, sn(Z,k') = a k) cn(Z,k')= — dn(Z, k^) к sn(z,k) 
sn(z,k)= — en(z,k) = Ue.) ipn 


L'équation de Lamé dans la nouvelle variable s'écrit : 


»"(z)- (h -v(v +1)dn° (Z, Ky) =0 


S y"(7)+ (v(v +1)=h+v(v « 1)? dn? (Z,k'))y(2) =0 


Dans la nouvelle variable, la fonction est de période réelle 2K' ou 4K' : 


p y(z +2iK(k')) = у(х) E ))= x) 
o Ж 
ou y(z+4iK(k'))Q у(2) lou y(+4KkK(k'))= y(z) 

Le procédé de construction des fonctions que je nomme Eic et Eis est donc finalement assez simple : 
- réaliser l'algorithme de construction d'un fonction périodique de période 2K' ou 4K' avec le 
paramétre k' 

- les coefficients trouvés sont les vecteurs propres de la matrice construite avec les coefficients k' 
mais de valeur propre égale à v(v+1}-h, et l'identification des fonctions Еіс et Eis se réalise comme 
suit : 


Eic?" [z + 2iK (e)? )= Eic?” (z+ 2ix (t), 2) Eic?" [z k2)= Ec?"(z7,k") 

Ес"! (z +4iK(k'), k? e Eic?” (z +4iK(k'), k? ES Ес"! GC Ec [Z tc 

But (z +4iK (cl k?)=E E (z + 4iK(k' ), 2) Eis?" (z, k?)= Es" (7 a 
2 


Sy 


Eis"? (z +2iK(k'), (k'),k?)= E SR (z + 2iK(k'), e) Eit (z, k 2 = ЕЁ," (z ; к?) 


S, 


Les propriétés de parité des fonctions de période imaginaire se déduisent des propriétés des 
fonctions initiales de période réelles, soit : 


z =K(k)+iK(k')-i ZR c -2K(k)-z 
Z => 2K(k ')- z z— 2K(k)—z 
Ec?" (2K(k')-2, k? )= Ec?" (2,k")=> Eic?"(2K(k)-z,k*)= Eic?"(z,k?) 


)- 
Ec" Gce)-24* ваа) о tic ene) а) Bebe 
Ex" Q(e)-z.*)- -Es Gt" )o Bit") s) Bn 
ESUS (2K(k')- 7, k")= NT же ‚к?)= Eis"? (2K(k k)- z,k?)= - Eis?" (z, k?) 


v 


Z, 
Z, 


La parité sur z entraîne les propriétés suivantes : 


K(k)+iK(k')+i7 =2K(k)+2iK(k')-z 

= 
Z>- | э 2K(k)+ 2iK(k')=z 
Ec?" (-7,k")= Ec?"(7,k")=> Eic?"(2K(k)+ 2iK(k')-z,k*)= Eic?"(z-2iK(k'),k?)= Eic?" (е) 
Ec?" (7, k" e Ec?" k )=> Eic2?""(2K(k)+ 2iK(k')-z,k)= Eic?""(z-2iK('),k*)= Eic?""(z + 2iK(k'),k*)=-Eic2""(z, k?) 
Es!" (Cz, )- E2 (Z, k? )> Eis?" (2K(k)+2iK(k')-z,k?)==-Eis?"*"(z - 2ik (k), k? )=-Eis?" (z+ 2iK(k'),k?)= Eis?" [z k?) 
Es?"(_7,k")=-_Es?"*(7, k" )=> Eis2"(2K(k)+ 2iK(k')-z,k*)=-_Eis?"2(z_2iK(k'), k?)=-Eis?"*2(z,k?) 
Construction des solutions périodiques d'équations modifiés de Lamé 


Finalement la construction précédentes nous a permit d'affirmer que les solutions périodiques de 
l'équation modifiée de Lamé : 


»'(2)-(h-v(v1an^(2.&))(2)- 0 


correspondaient aux fonctions périodiques construites sur le paramétre k' mais avec des valeurs 
propres égale à v(v*1)-h, h étant les valeurs propres obtenues avec le paramétre k'. 


Continuons ces transformation en considérant l'argument z-K(k)*iz', cela donne : 


= di A <i 


2= K(k)+ iZ = K(k)-7 


En conséquence modifions l'équation : 


qui posséde les mémes valeurs propres que l'équation ›'@)—(к-у(у +1їп!(@@,к))у@)=0/ avec des 


fonctions périodique construites avec : 


Ec". )= Ec?" (K(k')-z,k") EE, k?) Eak (ky Z ke) 
gs" Qa )- Bd k)-z,k?) вз") mu DT: 


) 


Si maintenant on revient à l'équation de Lamé originale г y"(2)+(h-v(v si? k))y(z)= 0: Par la 


transformation z _ K(k)- z il vient une équation modifiée de Lamé de méme valeur propre h: 


vafa -v(v «ie 22 EE 0 


ап? (z, k) 
Dont les solutions de période réelle sont construites comme suit : 


Ec?"(7,k2)= Ec?"(K(k)-z,k*) — Ec" (0,8) Ec?""(K(k)-z,k?) 
е) в) ageh se (t) ze 


Polynómes de Lamé, cas d'un développement fini 


Voici un tableau exprimant les huit types de développements possibles des polynômes de Lamé 
avec leur différente notation établie dans les divers articles fondateurs : 


Туре de Type de | Notation et Récurrence du | Notation et Récurrence du polynôme de | Normalisation NIST 
polynôme polynôme développement d'Arscott |polynôme de Lamé | développement de Fourier | Lamé selon NIST, formule dans 
de Lamé |de Lamé selon Arscott de NIST « A.Erdélyi.H.Bateman, HIGHER 
(Arscott) | (NIST) TRANSCENDENTAL FUNCTIONS VOL 
111, Lamé-Heun Functions, chapitre 
XV, section 15.5.1» 
s 2 pen 
1 1 uE" (2,2) ach, п>0 ||® PUE ijr) vc, A m0 e lons2) pour p=0 A +5 4, E 
У B. =4r + 1 2 P 
„ zn d r | po -(2n-2p-12n+2p+2)k? pow p21 p 
uE? (z, k? )= KEE (z.k) 7 --Qn-2r«2 2n «2r -1)€ uj. )- 4 4.84, Cos Deel 1 A, C 
ae Рі B,=4pl2-F) +54, > 0 
Normalisation 4 =1 2. UL 
1 = 
== „ (an- HIE NET ou A= 
Formule (18) avec v=2n 
2 2 a, =(2r+2)2r+3) d ls d e E 2 p-n 
Ea ser) mef n} n20 g Grille) sb nk uch, A 120 =у@л-2р-!дл+2р+4} Уды 
"^ (р ol y, = Qn -2r 2n «2r +1)? e zur, Dale A pour p=0 | "^ 
SE Jer Bn) sz e) aio) | 19, =2-"*2 E 
m E 2рн р d 4 (2 5 z À 0 
"NS и) рал) 8) pour pii оры 7 
Normalisation р =| А | 5-0 
| REES EES ou p 
Formule (20) avec v=2n+1 : 
3 3 n 2 a, = (2r +1)2r +2) " | d Ee i pei 
d | me (0,....n) п20 Ë per c nk nep.) 120 a, =5@n-2p-1)2n+2p+4)k Ya. =| 
геп P pea nl š с; p=0 
c ek j- T) As (n) y, = -Qn-2r* 22n + 2r +1)k p, в)-У [Ж n((2p+1)p) us 2-k = (2n+1Y2n+2)€ pow p=0 | 
m pa Driff) pour p21 Y (2p+1)B, 4 >0 
Normalisation As) l p=0 
y, = EE 2p+2\2n+2p 41) oug —1 
Formule (24) avec v=2n+1 i 
4 4 dE" А P dÉ ch - и) п>0 а, =(2r+1)2r +2) dE" É p me. d п>0 (2n+1(2n+2) pour р=0 
27 B, =4r Drift 207 e а„=\1 
e ' 5 с, š > (2n-2p+1Y2n+2p+2)K pour p2l|| 
dE”. (ы? )= an z DX " ,k) y, = -(n-2r* 2n 2r +1)k di "Ë H E аб, Je Econ Cosi (po) B, =4p"(2-K) 
Normalisation C, =1 y. =Lan-2p+2)20+2p+1) 
Formule (19) avec v=2n+1 
2 = (2r+1)2r +2 1 ï 
2 Z cd LE) тє[0,.„п-1) n»0 | ^ (rete ) al "kel sch, a) IO DEET DEET 
2 | B, =(2r +1Y(1+k?) ? 2 / 
$ =n- p-n-l -— 
cdE;, TUNI qz АЗИ" die) zl) y, Qn -2r on «2r +1}? cdi " (2, ‚|+ dn(z Y D, Sil (2p+1}) E pour p=0 DUAL ,20 
mi К Е 
Normalisation C, =1 @р+0-) pour р21 ои р, = 
1,7 npn tp je 
Formule (25) avec v=2n 
1 -(or л 1 p-n-l 
6 5 (6) тє[0,.„п-1) п>0 | |% Dana) . к IN vc, A) [Х| а, 7 -Qn-2p-3n 1p 4e ! x el 
T B,=(r+2} «Qr e Sich 2 — 227 
SE es Rat ўв SE | |, = а-о) зои) DOE 2p+2)p) | 1e Cr+2/0-F) ESO NN. 
А x 2p+2 
Normalisation B =1 y, 2 zn-2p «2p 1e p : 
с ou В, =1 
Formule (22) амес v=2n 2 
7 6 sdi y es n>0 | [e =(2èr+2)2r+3) sd ku me 0,.,n-1) п>) a, - (0л-2р-10л+2р+2) [ye [єс Cual 
K " B, =(2r +1} +(2r+2Y e 2 NECS 
КАЕ 2»: AE) | |, =on-20n+2r+ | Е [E )- a, Ў, al (+00) , a pur p=0 | Men 
m „= ` 
Normalisation p =] (р+к) pour р>1 ou Ç =1 
7,7 30n-2pns2p ye 
Formule (21) avec v=2n 
° Š scdE 4 LE] sch Al ii Ge wil fl `) пе), A) ml а, =(2n-2p-1Y2n+2p+4)C LEE: LES: al 
Е, В, = Qro) 2 Ga 
„ЇЇ Hl: dÉ » ы a u 2p+2)D,, > 0 
scdE] (56) "AM LU = п-п) «il " d di, ыў, ,Sin (272) b, (2p+2) e k ) А 
ои D, = 


Normalisation D,=1 


r, =>(n-2p)2n+2p+3J' 


Formule (23) avec v=2n+1 


Polynômes de Lamé, cas d'un développement fini à l'aide des polynômes de Tchebychev de première 
et deuxième espèce 


Les polynômes de Tchebychev de première espèce sont définis comme suit : 
Cos(pp)=1,(Cos(p))= Cos(pe)=T,(su(z,k)) 


et ceux de deuxième espèce comme suit : 


Blech =U „(Cos(p)) > Sin p* р) = Sin(p)U, (Cos(p)) = cn(z, kU S ED H 


Sin(p) 


Voici un tableau exprimant les huit types de développements possibles des polynómes de Lamé à 
l'aide des polynómes de Tchebychev de premiére et deuxiéme espéce : 


Type de polynóme de | Type de polynóme | Notation et développement de NIST par les polynómes de Tchebychev de premiére et 


Lamé (Arscott) de Lamé (NIST) deuxiéme espéce 
d 1 pen 
ir | ne| nj n2 uE” (z, k )= 4 +2. 4,7, I (sn(z, Di 
5 
2 2 ; pa 
sb, 2k | тєл n> see 2 LA, ipa (slz 


сЕ" kl welt, A vi Stolz DIS (sn[z,k)) 
^ 4 dE, y HI mel, n) n> az eras] Ee, b 
5 7 art. ne| NS n>0 о.) о), ЯЕ NE 
; i vk wech A1 n>0 EE "Latz AER (sn(z,k)) 
p 
d ° Dad |z, me! "a n>0 El) CA SE ab, onte 
8 8 n» 


sd? d пеј), A) nl sd dë ШЕ BIO A, Un sk) 


nd 


Equations intégrales des solutions périodiques 2K et 4K de l'équation de Lamé 


Les fonctions périodiques 2K et 4K vérifient des équations intégrales de Fredholn homogène dont 
le noyau se construit à l'aide des fonctions de Legendre de premiére espéce. Délibérément je ne vais 
pas prendre les noyaux qui sont présentés dans « NIST Handbook of Mathematical Functions, 
chapitre 29, Lamé-Functions », mais ceux évoqués dans l'article de 1940 de E.L.Ince, « V. The 
Periodic Lamé Functions » qui me paraissent plus commodes à manipuler. 


Fonctions de Lamé périodiques 2K : Ec/"(z) pour v£2n+1 


L'équation intégrale s'écrit : tt уу?) LA Pn tn teca) Vv#2n+1 


-2k(e) 


La constante se détermine de la manière suivante en se fixant à la valeur 0 de l'argument pour 
v£2n+1, ce qui donne : 


2m 2 
Vv#2n+1 asi ?)= а ) Comme Ec” [z k? )= dn(z, e > C,,Cos( 2) ф= 4742 am(z, k) 


P(0) | duE?"(u,k2) i 
-2k(e) 


вв) 


2K|k?) 
P (0) | du dn(u, Je Ee colpe) n HO Seb собно) 


-2K(k2) 


p=l 


Ec?"(0,k2) _ Ec?" (0,42) 
T C.P, (0) 


dE 


Ec?"(z, g)- dn(z, Je Ec, Cos( 2) pour v£2n«l 


p=1 


D'où l'équation intégrale : 


Ec 2" (z, g)- Ec.) [I sn(z, k)sn(u, k))Ec?"(u, k) 
T оу -2K 


L'équation intégrale peut aussi s'écrire de maniére équivalente en utilisant les propriétés de 
périodicité et d'anti-périodicité des fonctions de l'expression : 


Ec" z, k? )= EG") Гар (k sn(z,k)sn(u, к)Ес2"(и, Cl v=u+2K sn(v -2K, k) = -sn(v,k) 


TC, P. (0) Jk 
Comme Ec" (zi? )- Ee" Cz, г) еі Ec?" [y -2K DN Ec?"(v, k?) 
= Ec" [- z,k2)= Fe?" [z k2)= SM Jan) (k sn(z, k)sn(v Ec?" (v, k?) pour vÆ#2n+l 


On peut établir également que : 


Ec?"(z,k?)=d dn(z,k E 2рф ) pour v #2n+1 


Ec?" (z, k? ) = жые | duP, (k sn(z, k)sn(u, k))Ec?" (u, k) 
0 v -K 


Fonctions de Lamé périodiques 4K ` Ec"! (z) pour v£2n 
L'équation intégrale s'écrit à l'aide du méme noyau : 


) 
Ec?" (z, k? )= ail Jy Ta PU sn(z, k)sn(u, k)) Ec?" (u, k?) Vvz2n 


-2K(k2) 


La constante se détermine de la manière suivante par dérivation de la fonction et en se fixant à la 
valeur 0 de l'argument pour vZ2n, ce qui donne : 


2x(e) 
Ec" (0,42)= zm (e) n (0) In sn(u, k)Ec?" (u,k?) Comme P. '(0)= -(v +1)Р, (0) 


-ak(2) 


2т+1\ 2 
= 2 ок ) Comme Ec" (e )= dn(z Jc, on Er) = -amlz,k) 


Hais 
k(v+1)P.,(0) [au sn(u, k) E" (n, I?) "m 


-x(e) 


Vven Jr 


Ec" (0,42) 


k(v +1)P ©] Cos( SCH p+ у 
р=0 


Ec?""(0, pg 2) 
A+} СР, (0) 


=> sn(u,k)=Cos(p) et ami(e) 


Ce qui permet d'obtenir les trois équations intégrales équivalentes: 


Ec?" (z, EL dn(z Jc, cn бон) p=7-am(z,k) 


PAUSE 


) 
v Пт CP. a ( 
( 


Hai? 
Ec" (z,k?)= al ja PU sn(z, k)sn(u, k))E Ec? Quae) 
2т+1\ 2 4K|k 
TL j du bloc) vv #2 


2m+l1 (o. k 


2 en js 


Kle) 


Klk 


P. (k sn(z, k)sn(u, EIS Ec?" (uk?) 


Fonctions de Lamé périodiques 4K : Es "`! (z) pour v£2n 


L'équation intégrale s'écrit à l'aide du noyau suivant : 


e) 
Es?" (z, k? ) = А2" (е2) спб, k) le спи, ЮР, '(k sn(z, k)sn(u, k))Es?"*! (u, k) 


-ak(e) 


Et la constante multiplicative se détermine de la méme manière : 


2m+1 2 
д"ч(е)- = E. (o, k ) Comme Es?" (ie )= dn(z DX 2p+ Z ф= É -am(z,k) et cn(z,k)= Sin(@) 
P.'(0) n en(u,k)Es2"" (u, k?) à 
-2K(k?) 
ziel, — ? Es?" (0,42) _ L _ еы). 4 Se, 
pes z =e | yc 
p.'(0) | йи ote] n, sat P(0)[ de эңе) D, (ape il у yi 
-k(e) p=0 Ze p=0 


Ce qui permet d'obtenir les trois équations intégrales équivalentes : 


Es?" (z, k?) = an(z, D Y, se p+ w) 


p=0 


( SE Es (ok?) _ k? ick ofa ü yn ( ( ) ( ) 2 ( 2) 
Eelst Sy ? cn(z du cn(u,k)P,'(k sn(z, k )sn(u, k Es," и, k 
i z(v +1)D,P. (0) ,P. (0) -2K(k 
Es?" (o e) akle) 
E O Милл qu) 
v+ v+l 
2g?" (o, ie de) 
Es" 42 = Ler z, k) du спі u, k)P.'(k sn(z, k)sn(u,k Es?" (u, k? 
er) x (v -1)D,P...(0) p E )b (E sn(z,k)sn(u, E) Es?" (u, k?) 


Fonctions de Lamé périodiques 2K : Es" (z) pour v£2n+1 
L'équation intégrale s'écrit à l'aide du noyau suivant : 


к\ё) 
Es?" (z, )= 22 (en, A) In enu, k)P, (k sn(z, k)sn(u, k)) E. Es?" (u, k?) 


-ak(e) 


Et la constante multiplicative se détermine de la méme maniére par dérivation et fixation à 
l'argument 0 : 


2ке) 


Es?" (0, k?)= (ey > "(0) du en(u, k)sn(u, k) Es?" (u, k?) 
-2K(k2) 
2m+2( 1.2 Es?" (0,42) 2m+2 д ; 
À, (k )= зке) Сотте Es; Lu )= dn( (z, k nä (2p+2}p) 9-7 -am(z,k) et cn(z,k)= Sin(p) 
kP, "(0) fa en(u, k)sn(u, k)Es?"?(u 2) 
-2K(k2) 
EU 2. RN 2 
Zell, — Qu) - кр] 
kP."(0) | йи cn(u, k)sn(u, k)dn( (u DE pep), im (ob [ao Cos(p )Sin( EEN 
-2k(e) pm) 


DEIER 2E"? (0, Cl _ 2E52"2'(0,42) 


2E Es?" (0. i) 
kP,"(0)D,  xD,kP"(0) 


x D,kv(v +1)P (0) 


conme E)» vien (0) Are) 


Ce qui permet d'obtenir les trois équations intégrales équivalentes : 


Es?" (2, k dÉ dn(z ШР Р, зіп psal 


2m42, 2 2K DN 
ser Je QE) ot I" en(u, EIS (К sn(z, k)sn(u, k))Es Es?" í, k?) 
?) 


4K\k 
en(z,k) а cn[u,k)P.'(k sn(z, k)sn(u, k))Es2"*2(u,k2) 


0 


2Es2mt2 (0, k?) 
x D,kv(v +1)P (0) 
4Es2" (0,1?) җе) 


[4 m+ 2 
DWG DR] > 2; k) ш сти ‚ЮР, '(k sn(z,k)sn(u,k))Es? (u, k ) 


Es?" (2,2) = _ 


Es?™(z, k? ) -- 


Développement des solutions périodiques 2K et 4K de l'équation de Lamé par des fonctions 
associées de Legendre de première espèce et contructions des solutions de l'équation de Lamé de 
deuxième espèce par des fonctions associées de Legendre de deuxième espèce 


Les fonctions périodiques 2K et 4K peuvent se développer à l'aide d'une série de fonctions 
associées de Legendre de première espèce à partir des équations intégrales présenté 
précédemment. 


Avant de poursuivre ces calculs, j'introduis une formule bien pratique du théorème d'addition sur 
les fonctions associées de Legendre, donnée dans l'ouvrage de E.T.Whittaker et G.N.Watson < A 
course of Modern Analysis » paragraphe 15.71, et rappelée dans l'article de 1948, Erdelyi, « XXVII 
—Expansions of Lamé Functions into Series of Legendre Functions > (les deux formules diffèrent 
légérement mais sont équivalentes, et encore faut-il les modifier pour se placer dans le cas des 
fonctions de Ferrer, autrement appelées « On the cut » soit lorsque l'argument est dans l'intervalle 


[-1,1]). 


P, Ë: + J 7 i-e Cos(o)]= X иен] P (z)P” (t)Cos(mo) zie É 1,1] 


Z T(v+m+1) 


Étant donnée l'orthogonalité des fonctions sinusoïdales et le fait que les termes n et -n dans la 
sommation contribuent de façon égale à l'intégrale, il vient : 


ја Cos( (nw)P )P [z+ I PI coslo ))= [do Cos( (no)P. p [a+ 1-2 IP Cos(o))= эл T n1) ы) »(r) 


T(v n1) 


Soit sur la valeur t=0 : 


[ao ca no)P, je ` costo 1 bc no)P, je [1-2 Cos(o ))E2z rv- n+) ps (0)P"(z) 


"real | 


En dérivant la formule de sommation par rapport à ш, il vient : 


Sin(o Mi-z 1-P? p' (ar 4 Ji-z Wi- e Cos(o )- S пЦУ=т+1) ри); (њо) 


` " T(v ml) 


Soit l'intégrale définie suivante, puisque là encore les termes n et -n contribuent de facon égale : 


d'A [do silo) )Sin( (по)р, ра соо Jr Јао silo )Sin( (no)P, rh, IC? IC coslo )+ зат nt pr Р 


Soit sur la valeur t=0 : 


М2 [do Sin(o)Sin(no)P, rode L )=di-z j do si Sin(o)Sin(no)P, ri code, ele (0)P"(z) 


T(v +n+1 


Revenons à la construction du développement des fonctions de Lamé en fonctions associées de 
Legendre. 


Fonctions de Lamé périodiques 2K : Ec,” (z) pour v£2n+1 


En injectant le développement en série trigonométrique dans l'équation intégrale, et reconnaissant 
l'intégrale définie issue de la formule de sommation, il vient : 


Ec" (ue )= dn(z, desee, Cos( Gel et E pour v+2n+1 


PS 
sn(u,k)= Cos) k sn(z, k)-41- qn? (z, k) Posons w = dn(z, k) 


= Ec?"(z,k?)= <en | |, duP, (k (k sn(z, k)sn(u, k))an(u, k) E D» t S C, со 279) 


PS 


Ha jene Seeche 


> ке"(„)=^ SCH 
z 


| do eli Costo cost2po) = 2r E22) peser 


r 2p+1 
Comme + P 


Гаво) зя Ed 


А Ec?"(0,k2) [с а) aa; 
= вач) аа leide Al Se, ты) ier шы) 


Soit donc le développement suivant : 


Ec?"(z,k?)= о tr) et p= 5 _amz, k) pour vs2n«4l 


p=l 


вот) Y C) 


p=l 


Ec” oe) C, Г(у—2р+1 
Ec?"(z,k?)=2 КЕР (RAA) 8,71 p=0 6,9-0 p#0 v#2n et v#2n+1 


p. 


Ec?" (o к) " C, r(2n- 2p+ +1) 
E 2m k? =J 2n WV» p p? 0 |p?» di k 
Cn (z, ) C,P, (0) iar r(2n+2p +1) 2n ( ) 2n ( (г, ) 


р.0 


Fonctions de Lamé périodiques 4K ` Ec"! (z) pour v£2n 


En injectant le développement en série trigonométrique dans l'équation intégrale, et reconnaissant 
l'intégrale définie issue de la formule de sommation, il vient : 


pom (z, g)- dn(z, dër caen + H p= — -am(z,k) 
p=0 


aman 12) 2К(Ё) 
Edel du P, (k sn(z,k)sn(u, k)) Ec?" (и, 2) Vv #21 
k(v+1)r Cp (0) эе) 


Ec?" (o, k à axle) 
Kv zs СР, (0) At 
Posons w= dn(z,k) 


e Ec?""(z,k?)=- E Re Sc. Јар м? Соѕ(ф )kos((2p+ (el 


k(v +1)z СР, (0 


» al 


Ec" (z k?)=- 


ë Ec" (2,k2)=- 


du dn(u, Hlp (k sn(z, k)sn( u, k) D jc, nap] 


Comme а +1)p)=2 E Г(у -2р) pir" (o)p?e (w) 


Т(у +2р +2) 


т! 2| pem 
(o, k?) & * I(v-2p) P?” (0)P2P*\(dn(z,k)) 


Ec?" k? Ze 
"RE hh Saken 


De méme Ec?" (0,42)= pe +1X-1) Cu 


р=0 


Soit donc le développement suivant sur [O,2K(k)] : 


Ec?" (z, g)- ап (= DNE бә) 9-5 -amlz,k) 
Ес" (0,2) (0) р), 


p=0 


2m4li p=% 29 ) 
Ec (z 2\2 C, (o. ) r(v p p2rr!(9)p2P*!(dnl z. k у 2 2n+1 
C, íz, ) Hv DC P. (0) р (0 à To +2p+3) $ ( Р: ( n(z, ) vz2n et vz2n- 


EQ" (0,42) er. T(2n+1-2p) 
E 2т+1 k? d 2 2п+1 piri 0 |p2rn di k 
Ca (z, ) k(2n+2)C, P, ‚(0 NP AE TOn 1342 p) 2л+1 ( ) 2л+1 ( n(z, ) 


Pour ze [0,2K (6) 


Comme ce développement semble périodique 2K par sa dépendance à la fonction elliptique de 
Jacobi dn(z,k ), pour qu'il deviennent périodique 4K et qu'il possède les propriétés de parité et 
d'anti-périodicité adéquate, je le complète sur les autres intervalles par les formules périodiques et 
de parité : 

Ec?" (z, ?)- —Ec HE zk 2) pour 2<0 < parité négative 

Ec?" (2, k*)=-Ec deep, -ak(e ук?) pour еркі Jakle- parité négative par rapport à 2K ou anti — périodicité 2K 

Ec?" (2. k2)= Ес 2m (Mod(z,4K(k 2) k) pour z>4K(k?) S> Ea, a n 4K(k?),k*)= gc!" (oe) Vn e № c périodicité AK 


Fonctions de Lamé périodiques 4K : Ee "`! (z) pour v£2n 


En injectant le développement en série trigonométrique dans l'équation intégrale, et reconnaissant 
l'intégrale définie issue de la formule de sommation, il vient : 


Es", g)- dn( (z AED, 2p+ Z Q= > am(z, k) Posons w= dn(z,k) 


2т+1 2 axle) 
zéckt BUE). py ae (tnt ut) 


л(у +1)D,P (0) -ak(e ) 
— Es?" (o. k) 5 © "ао ) | | 
< Es, (z,k )= i =F DP Fame 0 zk JÈ Pare » Ge an u, k) en u, k)P. > '(k sn(z,k)sn(u,k))Sin((2 p +1)@) 


е Es" z, k2)= 


r ol “et a enk z,k) Ўр, ја Sin(p)Sin(2p (ele M Cos(p)) 


Comme 1% Sin(p)Sin((2 p мур) [i Cos(p))- Lee "n. Séi pir? (g)p2r" E 


E 2Es?"" (o, k) cn(z LE Tv- 2p) 
E 2m+1 k? Kn * 2 | 1 D. pr 0 Sa d k 
l WDE lte 200 + Dya c5," PT n) 


De plus Säck Co 


p=0 


On arrive au développement suivant sur IO 2K(k)] : 


Es?" (z, k ?)= dn(z,k In Б] 9-5 -amlz,k) 


p=0 
poo 


Ез?" (0,к2)= V C1) D 


2р+1 
р=0 
gen) 2E (0,2). en( 
Kv +1)D P, (0) sn 
2E (0,42) _ cn(z,k 
k(2n 2)D,P,, (0) sn(z,k 


(pen, PE OR ian) gemet 


DEI (2p+1)p, 1241-28) great Ai 


2p+1 Г(2п+3+2р) 2n+1 2n+1 


Je complète sur les autres intervalles par les formules périodiques et de parité : 


Es" k2)= Es Ge zk 2) pour z<0< parité positive 
Es2""(z2,k2)= —Es pna -2K(k?) k?) pour z€ bi ?)4K(k J< parité négative par rapport à 2K ou anti - périodicité 2K 
Es2""(2,k2)= Es 2m (Mod(z,4K(k? Luc) pour z> akle? je Ec" +n 4K(k2) k2)= Ec?" (ze) Vn eN < périodicité АК 


A première vue, il semble y avoir une singularité en z=0, toutefois l'étude de la limite suivante 
donne une valeur finie : 


_ vv +1)k E 
Lim Š спіс A P2P"dnlz, k))= 2 pur peu 
= k) 0 pour p>0 


Auquel cas on obtient une valeur de la fonction de Lamé en 0, cohérente avec le développement 


2m+1 2: 1 
se) 28 Ы о E ум Н sche s Ze 
Hv «TD, P..(0) ^ r(v«2) (v 1)P.,,(0) 


puisque ` p(6)=(v+1)P.,(0)- 


Fonctions de Lamé périodiques 2K : Es" (z) pour v£2n+1 


En injectant le développement en série trigonométrique dans l'équation intégrale, et reconnaissant 
l'intégrale définie issue de la formule de sommation, il vient : 


Es?" (z, k ?)= dn( (z Eo, s 22) Posons w= dn(z,k) 


2Es2"*2'(0 k ?) "ja 
Es” zk] =- k) |d k)P,'(k enz, El К))Еѕ2" (u, k? 
Sy (z, ) Steier, Ж cn(u, JA sn(z, )sn(u, ) Sy (u, ) 


2Es2"*2'(0,k2) 
z D,kv(v +1)P (0) 

2Es2"*2'(0,k2) 
x D,kv(v +1)Р, 


Comme jæ Sin(@)Sin((2 p + 2)p ИЕ ))= mnl Here 


o Es?" (2, k?)= - 


DS p-o 
cn(z, k) | йи dn(u, k)en(u, ЮР, '(k sn(z, k)sn(u, de D,,.aSin((2 p + d 


-2K(k?) p=0 


ә Es?" (2, k?)= - 


спі 2, k) ўр. 1% Sin( (e)Sin((2 p +2)p)P, p (iv coslo ) 


Jr 
AEs?"?^ (o, k) cn(z, k) 52% T(v -2p-1) 
E 2m+2 k? NS v E 2 2 D p? 0 p»? di k 
< BS, (z, ) DE v(v + DE, 0) sn Co 2! p+ ) 2p+2 Г(у+2р+3) " ( ) v ( n(z, )) 


De plus Es?" (0, k?)= SC 2p+2)D 


p=0 


2p+2 
On arrive au développement suivant qui d'origine possède la périodicité 4K : 


Es?" (2, k? )= 2" + d ф= 7 — am(z,k) 


p=0 


Es2"?(0,k°)= Y p, i1 "(2p«2 


p=0 


) 
ан cl 
d 


DE T(v -2p-1) „з 
Eelst À 2p+2)D p??? (9)p?»?? (qn(z k Vvz2n-«l 2n+2 
RE sis eh: = d co Жк л. 
AEs2"2 (o, k? ) en(z,k) Z T(2n41-2p) 
E: 2m+2 k? 2n+2 2 2 p2r+2 0 par d k 
52142 (z, ) D "(Ол +2)2л+3)Р,„.(0) sn(z, k) x p+ )D 2p+2 Tn t542p) 2n42 ( ) 21+2 ( n(z, ) 


A première vue, il semble également y avoir une singularité en z=0, toutefois l'étude de la limite 
suivante donne une valeur nulle: 


Lim aer ae) =0 Vp>0 


50 sn(z, 


Calcul des dérivées premiëres et secondes 
A partir des quatre développement établis précédemment , par dérivation on obtient : 


Eoee C, r(v-2p«1) 
Cp (0) 146,09 T(v*2p^1) 


р,0 
Ec?" (o, k ?) g 
k(v +1)C,P, 


v+l 


Ec?" (z,k?)= 2 P?P(0)P?2”(dn(z,k)) 


E $c, IP peque ais) 
) Z T(v+2p+2) 
= 


б 
ЕИ. 
gea P 


AE"? (0,42) e T(v -2p -1) „apat en(z k) „зә 
EUREN EGRE EX 2 2)D P?*(0 >K) p2p+2( q 
D, k*v(v +1)P.(0) = { РБ ы Г(у+2р+3) ” ( E "Weil 


rv- 2p) 2p+l cn(z,k) 2р+1 
+1)D. ET P” (dnlz, k 
) 2BH Ts T{v+2p+2) ' ( 1. z,k) v ( n(z, )) 


Es?" (z, k?)=- 


Les dérivées premières : 


ва) В Co Dear sy DM ap унааң. (etae st) 


Ec?" (0,42) poo 


p.0 


бө pog) "delt, Select к) (и ан 9?" (dn) 


Ecir(ak*)=-2 E P (0) < T(v -2p42) sn(z,k) 

Zeil L2 2Es;" M, k^) SS г(у 2 ) 2p+1 (v- 2 )cn? (z, Ю)Р »2P*\(dn(z, k))- dn(z,k 1+(v +Ven?*(z,k) \P2”*"(dn(z,k)) 
À a ) e XS Cr UD xp T(v MR j^ ol : ar ) ) 
see EE dër ap, Del sete) a) 


Et les dérivées secondes : 


(v +v +2)en?(z, k)dn*(z,k) k La (z ,KYdn?(z,k)+k'en?(z,k)))P2?(dn(z,k))- 
P*(0) -(v-2p+1)dn(z, + (у +3)cn? (z, Р, IP? (dn(z,k))+ 
+(v -2p* yv -2p+2}en (z, k)P?^ (an(z,k) 

Soen m (v +100 ect, Eltz, Ee sn?(z,kdn?(z,k)+k*cn?(z,k)))P2""(dn(z,k))- 
ee jer FEEL SD pit] ph + 

m dro p + (v - 2pYv -2p* lr? (z, РУ! (dn(z, DI 

geg d E (v + HI + 3)er (z, k an" (z, DÉI si (z, е + (v + Nen?(z,k))+(3v + DVD K)” (anle, k))- 
Es" (2, k2)= ER E (0) -(v-2p)dr(z,kJ8- ov +3) Jer (z, п) р, (анг, к))+ 
PUE AUS «(v -2pKv -2p en (2, k)P?^ (dn(z,)) 
(v +2)у + Alz, k an" (z, DÉI sn°(z, wet +(v + Der?(z, H))+(3v + 5)dn*(z,k)))P2(dn(2,k))- 
(2р+2)р,,, Dr (0) (v -2p-1)dn(z,k)3+(2v +3)en?(z,k))P22(dn(z,k))+ 
+(v-2p-1)v-2p)en*(z,k)P28(dn(z,k)) 


EQ, Р) p=% Cap Г(у-2р+ 1) 
GEI hui (z, DE 01+0,, Г(у+2р+]) 


Ec?"(2,k2)=2 


v 


AE (ok knlz k) ë 


E 2m+2 P = 
epe) Dale UE Wl: H 


p=0 


Récurrence sur les développements en fonctions associées de Legendre de premiëre espëce 
A partir des quatre développements : 


Ec" (oe) Cp T(v-2p«1) 
SCH £16 n 


p.0 


E 2 © 
2т+1 2 3 (o. k À < rv- 2p) 2р+1 2р+1 
Ес, [zk ) 2 {юш 9246 Tv«2p42)" (o); (dn(z,k)) 


vH е 


ckt Fire 


2Es?"" (o, k? ) $ E T(v-2p) Я cn(z k) 
Es 2т+1 k? 2 | 1 D. pir 0 , pir d „k 
5, (z, ) Hv + pP. (0) 2 P ) 2p+1 Г 42p42) v ( (22) v ( n(z ) 


4E52"*2 (0,42) e r(v-2p-1) cn(z,k) 
Es?" (z, k? 2p+2)D,,., P??? (0) —— Р? (an(z,k 
Rene "ës Оша)” O) 


La forme générique des développements est proposé : 


ckt ^r Dr ` (dn(z,k)) + "tie: 2p -1)P?" (an(z, Dai с, Ud A 
sere) Ў агза) x, „ e De, c, П, 
bi ue) MERS x, 70-2708" aie) ха р, D e 
sota) eS v р-у” цы) x, =p CREO p, Feste ү 
Pour la suite des calculs, j'utilise les expressions suivantes : 
— > еы PQ ET EIE 
Ee 022020) Rog tope 
E раан 


E(0-(v«)n.(0) — Eil zl +1), (0) 


Je rappelle les quatre développements trigonométriques : 


dy = 


< (v-2p)v+2p +1) pour р21 
B, =4p*(2-k2)-(2h-v(v 0e) 


e SEET HE 


Fonctions de | Type de développement Coefficient de la récurrence Formule de la récurrence 
Lamé 
périodiques 
mf, 12 Я Es 2 = 
Ес ) COTE Yc esto). v(v T Dk pour p= 0 Hs RE a.C, -0 
р=\ 


Y „Cap-2 +8,С,, *a Cy =0 


ee" |= (U) D, däs 
p0 


p, =(2p+2f(2-k?)-(2h-v(v+1)k?) 


y, Ch 2p-1\v+2p+2)k? 


a =h HBK (ET BEN 


n се) Ec" £)- ae.) С, coop 1p) а, =-(v-2p-1\v+2p+2)K" p "ect 
Pa 2 Y Cy +С +a Co Í, =0 
` 2-р EK -(2-v(v +1)E) pour p=0 
LE 
(2p+1f(2-2)-(26-v(v +e?) pour pèl 
?, =519 200+) 
Ез," (z,k?) Ex (e )e dalz, KE D, ip) a, = -(v-2p-1)\v+2p+2)K = ол 
ш 2 Y „Dap1 + ВО, +a D, 70 
ja 1 +)? а-у) Dour p-0 
È 
Ge h-r)-2r-vv +1) pour р>] 
?, == (v-2p)v+2p+l)' 
Es" (z, k?) 


BD, +a,D, =0 
Y Da, t DD t аР, =0 


Ce qui donne pour les quatre développements en fonctions de Legendre, une fois les substitutions 
faites dans la formule de récurrence : 


Type de développement Coefficient de la récurrence Formule de la récurrence 
ckt Te ree n) -(v+1v+2) pour p=0 BA +0, = 0 
ò = races E X, X. +6 X, =0 
P P P y-2p —(ү+2р+2у+2р+) pour pèl ” Le : S s. 
2 Ec?" (0,42) 
; E X, =2——— 
B, =4p2-K)-(2b-vlv+1)e2) T(v +1) 
v-2p+2 1 2 
E =- =--(v-2p+1Jv-2p+2)k 
un jr -2p+l\v-2p+2) 
вечае) yas rien tte) ¿= 293 EE L +Ó, -0 
i v-2p-1 2 pat BA pu +Ô Xa 13 70 
2+1) -Qb-vlv +e) pour p=0 Ec?" (0,42) 
B, = 2 X,22————— 
, XN | КГ(у+2) 
@р+0-]-(0-уу+)Р) pour p21 
v-2p+1 1 А 
=- =-(у-2рју-2р+1) 
Т QV -2p)v-2p+)) 
sorle) HEY x, rient) | 5 =Q 283 2035 _ 1142494243223142 | [Буй 8095-0 
SNZ, K ) ро у-2р-12р+3 2 2р+3 А „Xap1+B, Хы +8 Xa =0 
u ENEE EEN pour р=0 Rue Es?" "(0, 42) 
(2p+1P(2-k2)-(24-v(v +12) pour pèl кг +2) 
v- 2p+12p+1_ 2р+1,, 
=- ———=-—(у-2 2р+1 zi 
RIETSER d: SE 
вое) 2 IS т-р) (let) | $ = у+2р+42р+2_ 1 pili | [BX +044 =0 
бе 7 s) =-0, — v+2p+3]v+2p+4 "d 
"us IT P y-2p-22p+4 7 prr Fe EE 
H m+21 2 
B, =(2p+2f(2-K)-(26-v(v +1) uem 
2р 2p42 | l Í SE 
y- + + 
E=, pon =--(v-2p)v+2p+2)= R 
v+2p+2 2p 2 p 


Solutions périodiques 8K générales de l'équation de Lamé dites fonctions algébriques de Lamé 


Dans l'article d'E.L.Ince de 1940, < VII—Further investigations into the periodic lamé functions >, 
l'auteur décrit tout d'abord la récurrence générale pour un développement présentant la 


périodicité 8K et de parité positive : y(z)= УА (Э Ө} Еп injectant ce type de 
E 2 | 


développement dans l'équation de Lamé оп aboutit à Іа récurrence à trois termes suivantes : 


n=2h-v(v +1)? 


SEI +4r+3)2v -4r -Dk2A, „ n TEE SEI Ar -32v - 4r 5) A, = 0 


La récurrence à trois termes est équivalente à un systéme linéaire bilatéral dont la matrice est tri- 
diagonale. Là encore le développement est valide dés lors que h (ou n) est une valeur propre de 
cette matrice tri-diagonale. Une approximation du systéme linéaire consiste à tronquer des deux 
cotés la matrice à une taille donnée et de résoudre le probléme aux valeurs propres, en ordonnant 
ces derniéres de maniére croissante. Les coefficients sont les coordonnées du vecteur propres 
associés, dont on peut par exemple choisir comme normalisation que le terme central A, soit égal à 
À, 


n 22h-v(v 41) 


(4r 1f (2 € 4, 


- ; (2v + 4r - 32v - 4r UE A, + ; (2v + 4r -3JQv - 4r + K^ A, , = 47A, 


а,А + B,A, * y, A. = Qn -vv 1)? 4, 


rr) 


а, = - (2 +4r+3Y2v - 4r - D? 


B, =(4r-1f(2-k?) 
y, =- (Ov car 3o are s)e 


SA =[.. 4, .. 4, 4 А A, | M,A, =4A, 
= Au = ER Aa: A, А, А А, A, M, Aon = АЈА, 
By Q, 0 ca 0 | 
V na Bu D 0 
0 E i i 
0 y; Bi «у 0 . 
M,,. = 0 Yo Bo Œ 0 EC 
0 л B œ 0 
se 0 Ya Pa a, 
0 0 Jta В, Jd 
À valeur propre 2 
ргор AM. ep Ate +16 
А, vecteur propre 8 


A partir de cette construction on en déduit les deux types de développement paire et impaire, 
solution de méme valeur propre h et linéairement indépendante : 


r=+œ 


MORDI acod | 


г=—%0 


4r —1 
2 


оне) x G)- аз 


r=+œ 


F=—00 


( 


4r -1 
2 


onte) 


Il existe alors deux solutions spécifiques de période 8k pour lequel le développement bilatéral est 
simplifié, à savoir les valeurs telles que : 

2v =4р+1 

ou 

2v=4p-1 
Pour ses valeurs respectives, on peut montrer par la récurrence que le développement est limité 
dans un sens, mais il demeure illimité dans l'autre sens à savoir : 


2V=4p-1=4,=0 pour r> p= y(z)= X ace | A dE 


2v24p*1— 4,20 pour r> p= y(z)= Zacdl a d = X a,col( i не) 


r--p r=-_œ0 2 
Il est aisé de s'en convaincre à l'aide des relations de récurrence et le test numérique est plus 
délicat car les matrices sont limités en taille et le premier cas est facile à prouver à avec une 
matrice de taille limité mais pour le deuxième cas, il suffit d'inverser le sens des indices et de 
retrouver la propriété attendue. 


Pour autant cela ne nous indique pas en quoi ces deux valeurs spécifiques conduisent à un 
développement fini. Pour cela E.L.Ince développe une autre méthode en partant de l'expression 


analytique exacte de la solution pour „ = ! _ y(z)= Jan(zk)+ en(z.k)- 
2 


Dans cette hypothëse, Ince va donc proposer des développements de la forme : 


y (z) = dn(z, k) + cn(z, k) Y A, sn” (z, k)+ cn(z, k)dn(z, k) Y` B„sn” (z, k) 


r=0 r=0 


\ E J 


2v=4p+1 et рєм 


r=+œ T=: 


y (z) = Jan(z, k)— cn z,k ) > A, sn” (z, k)- cn(z, k)dn(z, k) > B sn” (z, k) 
r=0 


r=0 


— 


et 


y. (z) = dn(z, k)+ cn(z, k) dn(z, k)Y C n? (z, k)+ cn(z, k) D, sn” (z, k) 


r=0 r=0 


N — 


2v-4p-l et peN 


r=+00 r=+œ 


y (z)= Jdn(z,k)—cn(z,k) datz EIS C,sn" (z, k)- cn(z,k)> D,sn” (z,k) 


r=0 r=0 


0 


Et E.L.Ince montre dans ce cas que les développements sont en réalité fini, d'où leur nom de 
fonction algébrique. 


Solutions périodiques 8K spécifique de l'équation de Lamé dites fonctions algébriques de Lamé 


les deux algorithmes de récurrence pour la construction de deux types de fonctions de Lamé 
périodiques de période 8K. L'auteur recherche des solutions de l'équation de Lamé sous la forme : 


»(z) = Jen(z, k)+ dn(z, k)w(z) 
L'équation de Lamé pour la fonction w(z) s'écrit : 


»"(z)* (л -v(v +1)k°sn° (z,k))y(z) =0 


y) Jen raana u) = ut Gi DEE EIER 


sn(z,k) 


Posons H =h- er Lost), С (v I ac Е Jeste jete) 2 


sn(z,k) 


En transposant les deux formes choisies de développement de parité positive, comme suit : 


1 r=+co r=+0 
2v=4p+1 et peN-> y(z)- pe (z,k)= /an(z,k)+ cn(z, 2l y A sn” (z,k)+cn(z,k)dn(z,k) > B sn” (z, 0} 
Pty 


r=0 r=0 
et 


1 r=+œ0 r=+œ0 
2v-4p-l et peN »(z) = EC (z, k) = [dn( z,k ) + спі z,k ant k) > C sn?” (2, k) + cn(z, k) b D sn?" (z, Ө) 
e 


r=0 r=0 


Dans l'équation modifiée de Lamé, on constate que les développements sont finis lorsque p est un 
entier. Les deux développements donnent des fonctions paires de périodes 8K. Pour construire des 
fonctions impaires de méme valeur propre h, il suffit de remplacer cn(z,k) par -cn(z,k) avec les 
mémes coefficient, soit : 


1 r=+œ0 r=+0 
2v-4p4l et peN »(z) = ЕГА (2, k) = dan, k)- cn(z, 0 eh A, sn” (z, k)- cn(z, k)dn(z, k) > B, sn” (z, 0} 
Pty 


r=0 r=0 


et 


1 r=+00 r=+00 
2у=4р-1 et peNy(z)- Es. 2 (z,k)= Jdn(z,k)— cn(z, lat DY C,sn”(z,k)-cn(z,k) YD nik) 
pL 


2 r=0 r=0 


De plus on tiendra compte de l'équivalence des expressions suivantes qui permet justement de 
mieux voir la parité des fonctions mises en jeu : 


Per ӨКҮ Үү sr S (nn) + danse SE 


BE 1- kn (2.4) 
2 2 


On obtient les relations de récurrence croisée entre coefficients des développements comme suit : 
Ber pp s p reri 
4 
A,-B,-0 


Qr -2Y 4. 1p -2r(2r A1? д, -4(p-r Apr) 4, , Qr 2)B,, Qr (1 €? )8, 220€ B, , 20 


rl 


Qr -2Y B, lp - Qr tr «200 - € )8, -4{p-rXp+r+1)28,,-(2r+2)4,, =0 


r+1 r+l 


Développement fini Vr>p В, =0 et Vr>p+1 А, =0 


w=4p-1> p=? 
C;=D,=0 


(2r+2PC, a «lp - Qr +102 «2)8 -2r(2r+1)}C, -4(p° -r C, , -(2r+2)D,„ +(2r +1)D, =0 


ri rH 


Qr -2Y D, +ÍH -2r(2r «Y? - or +1)2r - 2), - (p? -2)€ D, , -(2r+2)C,. Qr 1)? C, 20 


rdi rz) 


Développement fini Vr>p C.=D,=0 


La finitude des développements fait que l'on peut les écrire dorénavant : 


r= r=p-1 
| ` Jan" (z, k)+ cn(z, k)an(z, k) Y Bn" (z, 0} 


2р+= e e 
< (rel Zu! i S 


2у =4р+1 et рєм 3 


z 
a sn Z.) x= r=p-l 
Es, 2 (z,k)= „241-2 2 p Asn” (z, k)— en(z, k)dn(z,k) Y B sn” c 
Wei r=0 


ЛЕ) 
2 


A КБ УГЕ) ce E 
m+— 2 2 r-p r=p 
Ес ? (z,k)= V2 К; > C,sn?' (z,k)+ cn(z,k) > Dan) 
SCHER = 


r=0 r=0 


et 


2v-4p-l et рєМ— 


Pour construire ces développements, je vais là encore me ramener à un problème de valeurs 
propres matricielles puisqu'il est facile de se rendre compte qu'une récurrence à trois termes 
correspond à la recherche des valeurs propres d'une matrice tri-diagonale. 


Si l'on sépare en deux blocs les coefficients A, et B,, ainsi que C, et D, alors il est facile de se 
convaincre que le problème matriciel consiste à construire une super-matrice 2x2 constituée pour 
chaque bloc d'une matrice tri-diagonale. 


On écrira symboliquement les deux systèmes matriciels : 


2v — 


2v =4p+1= p= 


M 


11 


matrice (p +1)x (p+1) 


matrice px (p + 1) S 


A 


A,2B,-0 


«А. + BA + Y; A, 4 + С.В 


rH 


* £,B, T ns = НА, 


Problème aux valeurs propres 


pea = pas HH 
4 

° G 

yi pi 

M, 0 .. 

matrice pxp | 0 
0 

0 а 

vi B 

M, Bi M 

matrice pxp | 0 
0 


Probléme aux valeurs propres | 


X -[c,.»,] 


г2 0 
C2B, t8; B, «n; B, , +a; A, + BA, * y; A, = HB, 
Développement fini Vr>p B =0 et Vr>p+1 А, = 
1 1 
- 0 de 0 
B al 0 44. 0 Eo бо 
1 1 1 n & e 0 0 
h B | 0 0 0 
0 ie x u M, 2 i S а e Ú 
0 7,2 Hz G 0 matrice (p+1)x p T» co NE 
Hl 0 y. Ў В! i a! Я fem 0 7] 5-2 € p-2 Ó p-2 
p- p- p- 1 1 
0 0 y p. | 0 ыз 0 nm, Epa 
о 0 0 0 
a 0 cm 0 0| e Ç. O 0 | 
Brody 0 0 0 m mn Xy O 0 
0 у, Bis со, 0 0 matrice рхр Qe Nr Ser 0 
0 у» Bp G> O » 0 mj Ep2 ns 
0 у Ba | 0 EE, AES 
Asa = [4 4, | 
Mi, M; Ара Ара 
= Н В, = E B, Í => Nombre de valeurs propres =2 p +1 
M, M» B, B, 
Xo pa T nl 
C. =D =0 
> 0 ds + BIC, ЕҢ СЫ + PD + єр, + тр, = HC, 
r> 
Ç; Б +e; D, +n; D, +a; C. + В2С, +У;С, = HD, 
Développement fini Vr>p C,-D,-0 
0 ... 0 | pg Zl Ò - à 0 | 
a 0 0 m & ¿& 0 0 
Ge га 77 M, 0 m 34 Ж 
Уз B3 О; 0 matrice рхр |.. 0 iau dw. 0 
0 yp Йб„„ Gy > Ò; Mrs Ve Gë 
0 Уа Bpa 0 0 np End | 
0 0 | FEE AE a: 0 
ar 0 0 m ë G 0 0 
TE s M, BEL s is 3⁄4 E 
qb Bis Qs 0 matrice pxp ` 0 туу uu бу; 0 
0 Уз» B; GE Ut 0 "ns 25 ёз 
0 Vi P L 0 0 т EN | 
calé 


C 
zl d D, = Io, o = Nombre de valeurs propres = 2p 


Pour la recherche des valeurs propres et des vecteurs propres, on ordonne par valeur croissante 
des valeurs propres et l'on normalise les vecteurs propres de telle manière que le premier 
coefficient du vecteur propre soit égal à 1. 


Valeurs spéciales des fonctions périodiques 


À des fins de comparaison avec d'autres implémentations, voici les valeurs spéciales en z=0 : 


1. mit 142 
2v-4p4l et pen ler 3 (0,4)=V2(4,+B,) = Es" 1'(0,&)= E (4, -B,) 
Pty Pt 


1 1 2 
2V=4p-1 et pen ln 2(0,k)=V2(C,+D,) = Es 2'(0,&)= Ex (A, - B,) 
Zë 2р-5 J2 


Pour le développement d'ordre général d'une fonction de période 8K, il vient : 


y*(z)= X aco [z dn 0) 


r=—co 


y (z- DIES dn veel! ai Jan 2) 


r=-œ 


2v =4p+1>4,=0 pour r> p= 


y (2)= Lace T dn 0) 


г=—0%0 


y (d= $X = "tard 1 dE 


2v =4р-1= 4,20 pour r> p= 


r=—co 


r=+p r=+p 


2»(0-X4 Oo > [£ 


2 


Tableau de quelques valeurs propres données par E.L.Ince dans l'article < Further Investigations 
into the Periodic Lame Functions » 


k=0.1 

(1) 2v=4p+1 m=0 т=1 m= m=. m= m= 

(2) 2v=4p-1 

у= р=0 (1) 0.275 

v=1+% p=1 (2) 0.42106079858305445 | 2.3289392014169454 

v=2+% p=1 (1) 0.6472914746853224 | 2.5983168305965796 | 6.379391694718099 

v=3+% p=2 (2) 0.9302172603145321 |2.9944478321576793 | 6.744804141786673 12.430530765741116 

у=4+% p=2 (1) 1.244338356787583 3.5278229121586446 | 7.223936253789767 12.897130559805731 | 20.48177191745825 

v=5+% p=3 (2) 1.5699630131322837 | 4.198309086712625 | 7.824381591160397 13.473901728975493 | 21.05039787094237 | 30.533046709076885 
k*=0.5 

(1) 2v=4p+1 m=0 m=1 m= m=. m= m= 

(2) 2v=4p-1 

v-A p=0 (1) 0.375 

v=1+% p=1 (2) 1.0089745962155616 | 2.741025403784439 

v=2+% p=1 (1) 1.729248688935409 4.375000000000002 | 7.020751311064587 

v=3+% p=2 (2) 2.4438189972283633 |6.459983068959071 | 9.290016931040922 13.306181002771634 

у=4+% p=2 (1) 3.153928595494932 8.637107257689303 | 12.375000000000018 | 16.11289274231069 21.596071404505036 

v=5+% p=3 (2) 3.8625805199677554 | 10.790225968477847 | 15.949363160436986 | 19.800636839563033 | 24.959774031522144 | 31.887419480032236 
k*=0.9 

(1) 2v=4p+1 m=0 m=1 m= m=. m= m= 

(2) 2v=4p-1 

v-A p=0 (1) 0.475 

v=1+% p=1 (2) 1.4210607985830541 |3.3289392014169454 

v=2+% p=1 (1) 2.3706083052819023 | 6.151683169403421  |8.102708525314677 

v=3+% p=2 (2) 3.3194692342588836 | 9.005195858213323 | 12.75555216784232 14.819782739685467 

у=4+% p=2 (1) 4.268228082541713 11.852869440194286 | 17.52606374621024 21.222177087841267 | 23.50566164321249 

v=5+% p=3 (2) 5.2169532909231044 | 14.699602129057682 | 22.27609827102451 27.925618408839664 | 31.551690913287405 | 34.18003698686772 


Exemples de valeurs propres et de fonctions périodiques de Lamé de période 8K (algébriques) 


Voici les premières fonctions de Lamé algébriques qui possède des valeurs caractéristiques h 
simples : 


1 
2у=4р+1 р=0 у : h jer) m =0 > Ec? JA k)= dn(z, k)+cn(z, k) Ens k)- dn(z,k)- cn(z,k) 


2 
= .Jan(z, k)+ cn(z, ln (z,k) js k^-41- k? + K" Jen Gil 
XE) en(z E) p DNE k?-41- Pe Jen n(z 
k) 


Bech, aped dye uh 100) I-k2+k* m=0> 


w BEES RS SS 


=a dn = n(z, 
2v-4p-l  p-l v : h te) k+k* m=1> 
Es? ке, k) )+[- k? + V1=k? +K" balz, k) 
2 


Les autres valeurs propres h ne possèdent pas de valeurs algébriques (racines d'un polynôme de 
degré supérieur). 


4) 
zc P ILE б, 0) 
) 


On peut aller plus dans l'expression pour la valeur particulière du paramètre К : 2 _ 1, pour lequel 
2 


E.L.Ince a également proposé une conjecture sur la valeur propre pour la valeur т=р, soit : 


pal 2у=4р+1 ы а т=р p- (4p+1)4p+3) 
2 2 8 
1 3 ; 1 
p=0 v = һ=——› Ес? = [dn(z,k)+cn(z, k) Es? = [dn(z,k)=cn(z, k) 
3 
Ec? = = ае) en(z, EI? -4sn?(z,k)-8 en(z,k)dn(z,k)) 
р=1 v=2 к= nu ` 
Es 2 =; dn(z, k)- en(z,k)(7-4sn*(z,k)+8 en(z, k)dn(z,k)) 
2 
2 
Ec? Е EC: )+ cn(z,k (13-36 sn°(z,k)+16sn“(z,k)-8 en(z,k)an(z,k)) 
E E 
Es2 - = ana, D ens —36 sn°(z,k)+16sn"(z,k)+ 8 сп(г,К)ап(,к)) 
2 
7 2 4 6 
ЕЗ = 2; Е Е) 347-1304 sn°(z,k)+1200 sn*(z, k)—320 sn*(z,k)- 
3 13 h 195 з 347 -16 cn(z, k)d dn(z,k)3 - 80 sn? (z,k)+ 40 sn*( z, k); 
p= у= = > 
2 8 7 2 6 
Es? he Det 347-1304 sn°(z,k)+1200 sn*(z, k)- 320 sn°(z,k)+ 
> 347 +16 cn(z, k)dn(z,k)p3 - 80 sn ?(z, k) k) +40 sn*(z,k)} 


Cette conjecture a été prouvée par A.Erdelyi, dans l'article de 1941 « XXVIII, On Algebraic Lamé 
Functions ». 


Dans ce cas de figure la récurrence s'écrit : 
Mere (4p+1)4p +3) P 


3 
2? sopla] 
| : : p(p+1) 


2v =4p+1= p= 


^a 


(2r+2) A, +{2p(p+1)-3r(2r+1)}4, -2(p -r41Xpr)A,. (-@г+2)В„, +(2r +158 .—rB,, =0 
(2r - 2Y B + +201) o Me 


r4l 

2(2r+2 B, +44 p(p+1)-3(2r «12r 4 2)B, -4(p-r{p+r+1)B.,-2(2r+2)4., = 0 
Développement fini Vr>p В, =0 et Vr>p+1 A4,-0 
2(2r+2 А, = 20r +2) B.,+{4p(p+1)-3(2r+1)2r+2){2r+2)8, -4(p-rXp r 12r 2)8, , 
22p(p+1)-3r(2r+1)}4, - 4(p- r41(p^r)A,, 
+2(2r +1)(2r +2)(21 +3)B,, 
+ (р Xr -2)-3Qr +1) (2r + DO 
-2[r«2(p-rXpr«1r^2) B, = 


h 2r +27 A. +2{2p(p+1)-3r(2r+1)}4, - 4(p -r -1Xp -r)A, , -2(2r +2)B, - 32r -1)8, -2rB,_ 


Mais je ne suis pas arrivé à la simplifier , afin d'avoir une expression simple des coefficients du 
développement. 


Fonctions périodiques de Lamé de période 4sK 


Dans ce cas la période est un multiple de 4K. E.L.Ince construit également les fonctions périodiques 
paires en utilisant le développement suivant : 


)- Y ACos((2r —0)am(z,k)) 0 = E p et s sont premiers entre — eux 


г=—0 


п z2h-v(v +1) 
Му-Ө+2г+2}у+ө-2г- De 4a n-Qr-oy oie )a +6 042r -lfv«0-2r- 2) A, = 0 


=0 


La récurrence à trois termes est équivalente à un système linéaire bilatéral dont la matrice est tri- 
diagonale. Là encore le développement est valide dès lors que h (ou n) est une valeur propre de 
cette matrice tri-diagonale. Une approximation du système linéaire consiste à tronquer des deux 
cotés la matrice à une taille donnée et de résoudre le problème aux valeurs propres, en ordonnant 
ces dernières de manière croissante. Les coefficients sont les coordonnées du vecteur propres 
associés, dont on peut par exemple choisir comme normalisation que le terme central A, soit égal à 
1: 


9 =P p et s sont premiers entre — eux n =2h -v(v + De? 
s 


rtl 


EE EE EA +(2r-0Y(2-k?)A, v 042r 1v -0 -2r « 2) A, , =n A, 


Sa, A, + B,A. * y, A 


r^ re 


=n4, 


r-l 


r 


a =-3(v-0+2r+2)v+0-2r-1K 


ФА, sl, A, … 4, Á À … А, ...] MLA, =nA, 

Se Au ER A A, 4, 4, d А, A,] M.A. = 7] Au 
Bno Q, O Dai the nnt vu. xe ` O | 
Y-n+1 Ps Q ny1 0 
0 


7 sioe vue 0 Yu D, a, 
0 a s du te 05 Wu. of 


À valeur propre A+ 4v(v + 1? 
d au & h = 
A. vecteur propre 2 


A partir de cette construction on en déduit les deux types de développement paire et impaire, 
solution de même valeur propre h et linéairement indépendante : 


DES Aco (2r -2 Jante) 5-9? азн (2, E 


r=-œ qu 


Il existe alors deux solutions spécifiques de période 4sk pour lequel le développement bilatéral est 


simplifié, à savoir les valeurs telles que: „Р = Qu  y-.P. 


S S 


Pour ses valeurs respectives, on peut montrer par la récurrence que le développement est 
unilatéral à savoir : 


vu > A=0 pour r<0— wt )= 4 coll Jante 0) 
S r=0 5 


у=-2_15 B =0 pour r<0= y'(z J- E acos (2-2 DD 


S r=0 


La récurrence sur les termes du premier développement pour „ - P est alors la suivante : 
S 


ET EE p et s sont premiers entre — eux п =2h -v(v De 
s 


—(r+120-2r-1 An «Qr -eP o - £?)4, -Qr 10 -r 1€ 4, , =n A, 


r+l 


— u +7,41 =n 4, 
=(r+1)(2r+1-20)k? 
= (2r-0ÿ(2-K°) 
=> - 
eA,-[4 A … … … … … 4 .] MA. =EN A, 
BRA ДИИ Les de a $c egeo] М.А ш = ALS 
(Де бу. O sam ts dee 10 
л A œ 0 
0 ... 
0 
M, = 
TES ftt 0 V n-1 Bui dy 
[0 … … … … … 0 7, В, | 
À valeur propre Маем, eo s EHE 
A, vecteur propre 2 


La récurrence sur les termes du deuxième développement pour ,, — 1 est alors la suivante : 


S 


v=0-1= Lo p et s sont premiers entre — eux n =2h -v(v +1}? 
s 


(2r+1-20)r 1)? B, Qr -0y (2-8 )B, - (C0 - r 12r -1)k?B, , =n B, 


Sa B. + В.В, +y, B,,-nB, 
а, =(2r +1-20)(r +1)k? 
( 


А 2r-0Ÿ(2-K°) 
=(r-1-0)(2r -1)k? 
к В, ] M.B, =7B 
e „а =[В В, В, 1 B,] M, nB,, 
B, o 0 0 | 
п В œ 0 
0 ... 
0 
M, = 
Ee, A Dr. Bo 053 
ID ue. gn S sau шй "ud Y, В, | 
4 1 2 
valeur propre | М. eh- À +4v(v +1)k 
В vecteur propre 2 


Les fonctions impaires sont linéairement indépendantes et construites ainsi : 


v=? => A=0 pour r<0= y ( z)- VA, зн (2r-2 Jon D 
s F s 


v=- P RB =0 pour r<0= y (z J- Ўзи (2r 2 DD 


M r=0 


Développement de Fróbenius des solutions de l'équation Lamé autour de la singularité 0 


Les solutions précédentes sont toutes utiles pour la résolution de problèmes aux limites et certaines 
correspondent à des fonctions propres d'un problème de Sturm-Liouville aux conditions aux limites 
homogènes sur un intervalle (celui de la période par exemple). En cela ces fonctions sont utiles pour 
la résolution des problèmes stationnaires de chaleur ou d'électrostatique, lorsque la géométrie du 
système de coordonnées conduit à des équations de Lamé, et la résolution s'effectue sur un 
intervalle comportant au moins deux singularités de l'équation de Lamé. 


Dans la plupart de ces problèmes, qui sont souvent à deux dimensions, par la méthode de 
séparation la fonction de la deuxième variable n'a besoin d'être défini que sur un intervalle 
restreint qui ne comporte qu'une seule singularité. Dans ce cas le formalisme des développement 
de Fróbenius est suffisant pour construire une solution. Et comme les fonctions de Lamé sont un cas 
particulier de l'équation de Heun, quoi de plus « simple » que d'utiliser les fonctions de Heun locale 
à chacune des singularités respectives (ces dernières sont en effet directement des développements 
de Fróbenius). 


Pour illustrer, voici les deux équations séparées du premier probléme électrostatique sur une iso- 
surface du systéme de coordonnées d'E.G.C.Poole de l'anneau aplati (Flat-Ring) : 


zi, на |» Jem (ua) (u) = 0 


du? 


ERIS AE G зне о)|өб)=о TC ma) 


Dans le problème électrostatique de l'iso-surface u= Ho porté à un potentiel donné, les bornes du 


domaine sont alors: “101 = et v e[- K(8. К(8)] pans le probléme électrostatique de 
l'isosurface у= vo porté à un potentiel donné, les bornes du domaine sont alors: 
uel0,2K(a)] е = vel K(B)-v.]o[ve KA) on voit donc que dans l'un ou l'autre 
probléme, il y a toujours un systéme de fonctions propres devant étre définie et réguliére aux deux 
bornes singuliéres de l'équation de Lamé, tandis que la fonction de la seconde variable ne doit étre 
réguliére que sur une seule borne singuliére. 


Dans le premier probléme le point singulier est O, et dans le second le point singulier de l'équation 
de Lamé est l'infini, car le point iK(6) correspond au point singulier infini lorsque l'équation est mise 
sous forme algébrique. 


La forme algébrique de l'équation de Lamé avec identification des paramètres spécifique de 
l'équation de Heun permet de construire < facilement > la solution de Fróbenius dans les deux cas : 


»(z)-0 


z 2-1 a?z-1 


d^y(z) 1| 1 | 1 \ a? dy(z) | 4h - a? (Am? -1)z 
dz? 2 dz 16z(z 1a? z - 1) 


Selon la notation des fonctions de Heun, notamment pour commencer, celle locale au point 
singulier z=0 : 


Fonction de Heun HeunG (a, q;&. B,y,8; z) 
por L. _ h g.l-2m ë HH] 
ei ag 4 4 


j y=ô=6=7>G+B+1=7+6+8=7 


Selon l'exposé très détaillé des 192 solutions de l'équation de Heun, par R.S.Maier en 2006, « The 
192 Solutions of the Heun Equation », pages 15 et 16, on peut établir les deux solutions locales 
régulières en 0, qui donne donc la solution de l'équation de Lamé sur l'intervalle u €[0,u, ] lorsque 
uo< K(a) : 


THW, fa G i Jem oso) (= 0 avec uel0,u,] 


h 1-2m 2m+1 1 1 


H 


4a? ” 4 4 '27/2 


H (1) нато A , EEN => H (0) = 0 ët H;'(0)=0 


І+а?-ћ 3-2m 2m+3 3 1 


° H 


Да? 4 4 272 


EEN а) 1 (0)-0 et H,'(0)> 0 


Dans d'autres problèmes aux limites extérieurs, u est dans l'intervalle u€[u 2K(a)], si uo >K(a) alors 
l'expression du développement est le même du fait de la parité de la fonction elliptique de Jacobi 
par rapport au point u-K(a) : 


CH. б азила) нео avc uel[u.2X(a)] 


1 h 1-2т 2m+1 1 1 
HIE n m " egre et H,(0)-0 


1 1+a?-h 3-2m 2m+3 31 , | 
А зп (u,a) 


sn*(u,œ)= зт (2K(a)-u,a)= 


=> H,(0)=0 ег Н, (0) 0 


H H $us 


EE xj 7 4 723 


Les récurrences de ces deux développements autour soit de O soit de 2 K(œ) s'écrivent alors comme 
suit lorsque l'on se réduit uniquement au cas m=0 : (elles sont plus efficaces que d'avoir recours 
aux fonctions de Heun implémentées sous Mathematica) : 


H(0)#0 et H,(0)-0 
ME 1111 RE a 

HE = gas) 2 2 (u,a) 
Жей йы: E = An - Av (i-o A, + (3-4) a" 4, , 

2 8(n+1)(2n+1) 
H,(0)=0 et H,'(0)#0 

l lza2-p 3331 

није nne A 23 iiber nta änt (u.a) 
yep йе КУ Жин”) n]B,-(1-4n) 028, , 

6 Sin +1)2n +3) 


En revanche pour les deux cas u€[0,uo ] lorsque uo» K(a), ou u€[uo 2K(a)] si uo «K(a), pour l'instant 
je n'ai pas trouvé de développement de Fróbenius convergeant car l'intervalle contient au moins 
deux points singuliers de l'équation de Lamé. Par contre dans ce qui suit, je propose une 
construction approchée d'aprés des fonctions sinusoidales hyperboliques. 


Construction des solutions de l'équation de Lamé sur un intervalle [0,xo] avec xo » K(a), monotone 
croissante ou décroissante 


Si l'on résume le probléme de la fonction de la deuxiéme variable, cela revient à trouver la solution 
de l'équation de Lamé : 


y"(x)+ (h -v(v + Ja?sn°(x,a)}y(x) =0 


xel0,x,] avec x,>K(a) et y(x)20 et be 2, SES 


ou le problème équivalent : 


y"(x)+ (n -v(v4 Da*sn*(x,œ))y(x) =0 


x e [x,,2K(a)| avec x, < K(a) et y'(x)<0 2 ОДЕ 5i om 


Dans le cas du premier probléme lorsque xo est inférieur à K(a) ou pour le deuxième problème si хо 
est inférieur à K(a), on a vu que les récurrences des deux développements de Fróbenius 
conviennent autour soit de 0 soit de 2 K(a). Il s'agit des développements des fonctions de Heun. Il 
s'écrivent alors comme suit : 


H(0)0 et H,(0-0 et H(x,)=1 


l h vvHll , NO 42034 
Hemo 3. d 3^ 9 jg (sa) LAm (x,a) 
Hi(x)= = 


à 1 h vv#ll TEN T 

алб, ы? 7 тта) 2,4 (x,a) 

EEE RBS AM. +(v(v +1)-2(2n-1)n-1)a24,., 
2 2(n+1)2n+1) 


H,(0)=0 et H,'(0)#0 et H,(x)-1 


n=+0 


2 
vien. E 1 s tel sn(x,a) Y. B,sn" (x,a) 


dei 2/2 22 
Biet: 1 1+a*-h I-v v+2 3 1 ` CM 
банана, E ZE y emsa) vna), Ba (хь) 
ке ждей p, lante) JB, «(v(v--1)- 2n leg. 
6 2n+1)2n+3) 


Mais ces développements pour les deux cas џє[0,џо ] lorsque uo» K(a), ou u€[u.,2K(a)] si uo «K(a), 
ne sont pas valables, car l'intervalle comprend une singularité de l'équation de Lamé. 


Ici, je propose une méthode de résolution approchée des deux problémes, sachant que le 
paramètre v est toujours supérieur ou égal à -1/2 et que les valeurs propres h sont toutes négatives 


et rapidement de grande valeur, par exemple dans le cas de la résolution du probléme de Sturm- 
Liouville pour la détermination du potentiel électrostatique sur l'iso-surface u=Cste du système 


de coordonnées « Flat-ring d'E.G.C.Poole » 


De plus les valeurs propres h sont croissantes en valeurs et deviennent vite plus importante en 
valeur que le terme jacobien. Dans ce cas on voit vite que l'équation de lamé est alors trés proche 
de l'équation d'un sinus hyperbolique, forme effectivement obtenue par la résolution numérique 
des problémes différentiels : 


h--A«0 | A»»v(v4lasn'(x,a) y" (x)+(h-v(v + ови (о) (к) 20 y"(x)-4 у(х)=0 
> у(х) ACosh| [7 tv(v + Des (x, a)}+ вета +v(v + Da sn (s,.a)) 
»0)=0 ` Mal y'(0)=0 EE 


4A Sinh Aevi Nm) ou ҖЫ Coshlx favi Her (s.a) 
Sins A +v(v + от (ха) Coshlx, A -v(v + Das” (а) 


Dans ces conditions proposons la forme suivante des solutions des deux problëmes : 


yG)-(ev(vete se sa)()20 лија (х,а) G 
SinhlxV Ï x 


x0)=0 7'(0)#0 ` all ` ek Sab Sid »(0)=0 x(0)#0 y(x,)=1 


_ Sinh fix _ [F Cosh mi _ Cosh wi ic 
deg Sinh уй d 5 I Д Sinh x VA Е zm Cosh уй Б 
el {слон Ta) LE rl JI 
y'()= И) а) cet) ` eg ` gll Pede) 20698)» sg") 
э yG)- [oet озо (а) а) Pils 2r Ge) fe") vl DG) Gs а) Ае) 
= (ут) X Ge) P Gde G9 re")? GG n) - Gs roe) 


= f(x)g"(x)+2 f'(x)g '(x )- v(v +1)a? (sn? (x, a)- sn cke a)f 0206 g"(s)+2LC fr )-v(v+1)a*(sn*(x,a)-sn*(x„a))g(x) = 0 


JF: У(0)= f'(0)g(0)+ f(0)e (0) 0 g(0)+0 Zi »(0)= f'(0)g(0)+ f(0)e(0)- 0 g'(0)=0 
g'G)e 2i Cotanil pe (x)-v(v +1)a? m (x,a)- sn (x, a)je(x)- 0 x) NI: Tail pen (x)-v( +1)a? m (x,a)- sn (x. a)je(x)- 0 


Je propose d'utiliser maintenant une méthode de perturbation de l'équation différentielle sur la 


2 
fonction g(x) dont le paramètre de perturbation € serait ç= v(v+1)a* | 


Soit l'équation différentielle suivante du problème aux solutions recherchées impaires : 


Ge 
= e eene sl) Nu" а E ,)=1> к) о 
8:(%)= 0 


Prenons la première équation différentielle, nous avons pour la première fonction go(x): 


g "(x )+2 JÑ Cotanh| Zb (х)=0= И Log(g,'(x))- een 22 пови i.) 


Log(g,'(x))- C, оз) g, (x)- US gx) C, C,Cotanl| ZS g (x)=1 


Pour la seconde fonction, l'équation différentielle est la suivante : 


g (x)=1— g,"(x)+ 2 À соат Zb, 2\/й (sn (x,a)- sn’ (x,a) 


>g (x) = Á+ BCotanh d х |+ solution particulière 


Une solution particulière doit être trouvée, mais il se trouve qu'il y a une relative difficulté dans la 
recherche d'une intégrale particuliére. Plutót que de vainement chercher, j'ai tenté finalement avec 
un certain succés de développer le second membre autour de la valeur x=x0, cela donne une 
équation différentielle comme suit : 


g," (x) WF Cotanh| х, (х) = АЙ (х-х, )зп(х,„@ )еп(х,„а Malz) 
La solution générale de cette équation est la suivante : 


EE ee M2) 2x A –4хх,4 окт) V7) 


Comme le terme en cotangente hyperbolique doit rester fini, il vient la valeur de B : 


B- snx, a )en(x,,a)dn(x,,œ) 


PNA 


Et la fonction doit s'annuler en x=x0 : 


dig mi aer onis a) | | сол, 


Уй 


D'où la solution de 9:(х) : 


g.(x)= — kn, a) Ë: Xo J X; Nï CT, SL i} xy Ji (сокту Xo |: Сап} x) 


Soit l'équation différentielle suivante du problème aux solutions recherchées paires : 


gx)-1  g(0-0  g"(x)+ 2 Tan xe) - A m (x, œ)-sn?(x,œ))e(x)=0 


'(0)=0 
mO e dtreg()=1>a(s)=0 


g(x)= g,(x)+ ë g(x)+ 8° g,(x)- = T 


Le traitement de la solution est relativement identique et nous avons d'abord go(x)-1, puis pour 
91(х) en développant le deuxième membre de l'équation différente autour de хе, il vient : 


g," (x) Anl x e (х) = AA (x x, )sn(x,, a n(x a Jin(s. a). 


La solution générale de cette équation est la suivante : 


g,(x)- А - Brands miscens rts.) l 2x! - 4xx À и? zwi] 


Comme la dérivée première de la fonction g1(x) doit s'annuler, il vient la méme valeur de B que 
précédemment : 
B- snx, a )en(x,,a)dn(x,,œ) 


PNA 


Et la fonction doit également s'annuler en х=хо, ce qui permet d'obtenir la valeur de А: 


d vis aene an. f тал EE re 


D'où la solution de 9:(х) : 


a (x)= sn(xo, a jenes. cz) (rann ДЬ таяй, ) xà) Ji X, fe x: WE Tanl ix) 1] 


À 


En résumé on aura donc la construction des solutions suivantes, sachant que le développement 
approché perturbatif est valable quelque soit la valeur de x, dans l'intervalle [0,2K(a)]. 


x e [0,x,] et 


<` 


x e [0,x,] e 


x e [0,x,] et 


x e [0,x,] et 


X, < K(a)> 


X, < K(a) 


x, < K(a)> 


X « K(a) 


H(0)+0 et Н,(0)=0 et Hill 


h +111 re 
gel $ E 3 T ER sn'(x, d ZA (xa) 


EE мол si (x,,a IS Уло х,а 


H(x)- 


gi ag 2222 


h 4n (1+0 ]- hJA, +(v(v -1)- 20n-T)n - Uer A, , 
а Aaa 77, Ara) dari ' К 


2 
Hat e оз ы Ез күт 


Ni 
g(x)= эда} T аја Had! ei — x] x) Ji ss x, ME em 


1 


MT fa (v +1)a?sn*( Tad) 
= Cosi Ja +v (v +1a*sn*(x,,a | 


-e g (x 


H,(0)=0 et H,'(0)#0 et H,(x)=1 


2 п=+0 
ЕЕ 2 m < k иа) sn(x,œ) У Bsn” (x,a) 


H,(x)= 4a? 2 2 2 2 - pe 
i [a^ ship v42 3 1 ; IN 
sn(x,,a )HeunG| - pre A 9515 pp (х,а) sn(x,,a) Y` „sn (х,а) 
п=0 
l-h+a? (ол-1) (1+) DEET GEN 
B,=1 B, = B, nd — 
6 2(n+1)2n+3) 


2 
H,(0)=0 et H,(0)20 e H,(x)=1 Urb i-Asv(vele (ха) h=-4<0 


ni 
g(x)= sus, entro Hd) ls X; е X М, Coin n 1} x, J kal? x ү Cota "| 


À 
ч Sinija eer] TUNE 


Sinhlx, үй +v(v + Da" sn (xa) 


Pour les problèmes sur l'intervalle [xo, 2K(a)], il vient pour les solutions approchées : 


DE er Bra а AH) e ME snn ш 


ni 
хє[х,.,2К(а)] et lou > ss) HP DS, "el ix f xA )+ Vi (x Ah МТ ixi] 


x, « K(a) (ок( 
.. Cosh(2K(a) )-3 1v (v +1)a sn (x, 0 ) | ЕЕ 
` Tosiek) =x) Vill +v (v +1)a? sn (x, œ JY | к ER | ) 
A H2K(a))=0 e H,'QK(a))#0 e н.) te EE ee 
x, 2 Klo 
ref 2K(a)] et {ou , i aia) оь ss „Мснин |) (отд, 3x )- com iJ 
x < K 


NEL OE ни а) a) экы) 


Emm x) Hi (v Tos (х,а )) 


Tandis que pour les solutions exactes, nous avons : 


H(2K(a)s0 et H'QK(a)-0 et Hj(x)-1 


1 h V vil 1 1 2 NS 2n 
H ; 4 
vol 5. a 32 yy visa) _ > sn" (x,a) 


TE h : V QE L n (x, a ) YA sn?" ( (xa 


xe|[x,2K(a ) et x, > K(a) > H(x)- 


a 433! 2! 2 7272 


aap gaba q -l)a + (t )-2 25 1Yn-D)a24,., 
l m EH 2n+1)2n+1) 


H,QK(a)-0 et H,'(2K(a))#0 et H,(x,)=1 


2 
were! L = 7. =k E) sn[x,a) Y` B,sn” (x,a) 


xe|[x,2K(a ) et x,» K(a) > H, (x)= 


п=+0 


sn(x,,a LP Eng ск” рне) sn(x,, a) у Bsn” (xpa) 
п=0 
1-л+а? [na (ie? )-4)8, (v 1)- 2022 ов, , 

‚к жы 5; 2n+1Y2n+3) 


чө 


Autres développements de Fróbenius des solutions de l'équation Lamé autour de la singularité 0 


I| se trouve que certains problèmes aux limites présentent également la nécessité de 
développement autour du point singulier O pour l'équation différentielles suivantes : 


d’@(v) l | (m je "Zei ve[ vovo] vo <K(B) 


dv? 4 cn? (v, B) 


>V = iv < 277 G 1 Jersm"(7.a)j6()=0 AMES Ga) 


Dans ce cas au premier abord j'utiliserais les deux développements précédents en substituant la 


2( > 
valeur 5" (v.a), Cela donnerait : 


h 1-2m 2m+1 1 1. sn/(v,B) 


oi 4 ' 4 '2'2 en (v, f) 


6e неш] i | > @,(0)#0 er @'(0)=0 
a 


1 l1+a*=h 3-2m 2m4 1 sn 
cn(v, В) a 4a 4 4 22 cn(v,f) 

Mais il se trouve que la récurrence de ces deux développements conduit pour certaines valeurs de 

vo à une divergence nette de la série car le terme sv (v.P) est tout simplement supérieur à 1. 


cn (v, В) 


Dans ce cas il convient d'utiliser d'autres développements de Heun équivalents au point singulier O. 
C'est encore R.S.Maier qui nous sauve dans sa publication de 2006, < The 192 Solutions of the 
Heun Equation », pages 15 et 16, avec les deux développements suivants : 


eo (1-3) * tnc CUT 841,y,d - B4 m 


=] d= 


өд) x"(1-3)*" Hana] ei 4+(у IE (v Jes y 842) y)a . ç PE ecd d 
а- a- x- 


Ce qui donne avec la substitution des paramétres de Heun correspondant à l'équation de Lamé : 


Т 1 h ï 1-2m j 2m+1 $ =g 1 a 1 
= E = H = 9 7 
ai а 4 x 2 a-1 B 
2m-1 П 
ө(х)= (1-х) + HeunG, ы md ш т; X 
p sp que usd 
1 3-2m П П 2 
@ (x)= (01-х) i HeunG, t š ente 2 773 mo m; X 
p 88 4 ? 4 2 x-1 
Une fois injecté la valeur de x = sn*(7,œ)= sn (v, B) їс 1 X = sn? (v, В), il vient 
cen? (v, B) cn'(v,B) x-1 


à une constante imaginaire prés : 


1-2т ï 
Gen (нед) 7, 2552. E Z i vs 
2н 1 7+2h-6m+2a° 5-2m 3-2m 3 
@,(v)= зп(у, B en 2 неј zr 71 a í m, = ms) 


Les récurrences de ces deux développements autour de 0 s'écrivent alors comme suit lorsque l'on 
se réduit uniquement au cas m=0 : (elles sont plus efficaces que d'avoir recours aux fonctions de 
Heun implémentées sous Mathematica) : 


@(0)z0 et @'(0)=0 


1-2т 


okt? 1 1+2h-2m 1-2m 3-2m 1 25 


напо. ТЕ : "cwn ips m;sn'( X Ast 
_1+2h 


iaa 4, Mehmet p) (5-161416 )?4, , 
4 
0 


8(n+1)2n+1) 
6,(0-0 et Ө, (0) 


1-2m i 1+ a 5 40 
0,(v)- sn(v, Ben Í = E i L 2s S Ki mhi m, sn Чу, ДЕ У Bn É E sn” 
В 88 4 A" 
+ s e а 27+2h-2a?+8n*(1+ В°)+ 4n3+ 8° )4, +(1-16n?)824, , 
12 8(n+1)(2n+3) 


Etant donné que Y° < K(B), il n'y a pas d'autres points singuliers de l'équation dans l'intervalle 
d'étude et donc pas de nécessité de faire appel à une technique perturbative pour trouver une 
solution approchée. Comme les deux solutions sont soit paire soit impaire, elles s'applique autant 
sur l'intervalle [0, vo ] que sur l'intervalle [-v,, vo J. 


Pour les problèmes suivant les solutions découlent immédiatement : 


d’@(v) [ | б 3 s ten) et)=o v elv] v,<X(8) 


dv? 4) cnm(v,f) 
]-2m 
сп 2 (v, B)HeunGi J. : m m 2". 2 m Ll m; sn? (v, D) 
9-1 e(0)-0» 69 )- — I 1+2h-2m 1-2m 3-2m | 
2 H ‚1 
сп ? (v) «тб э 857. Aa ou» m;sn m) 


zs 1 7+2h-6m+2a° 5-2m 3-2m 3 
2 HeunG 1 
sn(v,Bjen 2 (v, В)Неип EN 3: ; оса 0» m;sn°(v, 2 
1-2т 


S 1 7+2h-6m+2a° 5-2m 3-2m 3 
sn(v,, Ben 2 (v,,B ER sp bu doa x 


6(v)-1 @(0)=0— e,(v)- 


de. 2 


Construction d'une première catégorie de solution locale régulière à l'infini par un 


développement de Fróbenius 


La solution locale régulière à l'infini, de l'équation : 
2 2 
d av) 33 @ ЗЯ sn (v.f) 6(v)- 0 
dv 4 cn (v. B) 


F =iv Pop), [, б Lescht Ge ciel: ERO 


s'obtient également par un examen attentif des diverses solutions 1 et 2 donnée par RS Moler en 
2006, « The 192 Solutions of the Heun Equation » en pages 24, 25 et 26. 


Remarque : dans ce qui suit la variable x de l'équation algébrique figure l'expression : 


sn*(v. B) _, 1_ cn°(v, В). z 1 cn? (v, В) 


cn*(v.,B) x sn°(v, B) < 1—х ` en (v, B)+ sn? (v T" 


x =sn"(V,œ)= 


On peut donc mettre à profit la transformation homographique 1/(1-x) qui figure justement parmi 
celles donnant les fonctions de Heun équivalentes, pour obtenir une expression « simplifiée » de la 
solution de Fróbenius. Mais revenons à l'expression attendue des solutions régulières sur 


2 

l'intervalle : e vo K(8)]=> xe | œ, en Les 4 expressions des solutions 1 et 2 au point 
Vo B 

singulier infini données par R.S.Maier sont les suivantes, soit régulières soit non régulières à l'infini. 

On remarque que les solutions #1 et #2 sont interchangeables par la substitution : a<->B : 


Fonction de Heun point singulier z =œ 


WEE arda se a 8j sg raa- Past = 
z 


a a 


Solution #1 — 


| l q+@la(& Y ó +1) a +ó Jj. f 
1-0. а—1 1— 


E SÉ :B,B ét geet 
a a z 
Solution #2 — 


а—1 


ËM а+ Bl] -y -51)- B «o os ё+1,Й-— air J 
Em =Z 


Or nous recherchons une solution régulière. Dans le cas m=0, on a : 


N | = 


Alors les solutions #1 et #2 sont parfaitement identiques et toutes deux réguliëres à l'infini, ce qui 
nous amène à proposer la forme : 


v=iv > ZS) 4 d ü Je esaet )-o deb uero up) 


cn" (v, B) 
_ dek) [nef тезе) = Me. velo kif) 


eost nte) 


ort) ten | a Т 2: E 5 ‚сп? (v, 2) 


— 


<0 et 


Lorsque m> 0, alors : a = 1-2m B- 2m+1 > 0. Et dans ce cas la seule solution régulière à 
4 run 


l'infini est la solution #2, soit pour l'équation différentielle suivante : 


d’@(v) l | |» re eO avec velv,,.X(B) 


dv? ?(v, B) 
OR [men gd" E , (2m+1ÿ(1+a°)- Ah 2m+1 2т+3 |, 1, atA) 
Solution #2 — v.g) 16 4 4 2° sn'(v.p) 
2 2 2 
Lees = a? (2m+1V1-3a*+2mB*)-4h 2m+1 2m+3 у 1. , 
(cn 6,8) : 4 eG[ pU 160? бо а Ле (v. B) 


La récurrence du développement autour de l'infini s'écrit alors comme suit lorsque l'on se réduit 
uniquement au cas m=0 : (toujours plus efficace que d'avoir recours aux fonctions de Heun 
implémentées sous Mathematica) : 


d'e(v) [ 2 "äert avec v ch, K(8)] ауес @(K(B))=0 


dv? 4 cn'(v,) 
2 n=+% 
ov) = entem | а? 1-4h-3a* 1311. PP) cts 
p" deg^ ‘44 2 Z 
3a?+4h-1 (3a2+4h-1+8n(2n+1)2a?* —1))4, 6 -16n 16 fi- a 4, , 
А 21 A, =—— 4 Аа = 2 2 
16a 16а (n+1) 


2 WW 2 
HeunG, - 1 4h ai 1 3. I. cn" (v, p) 
cn v,p p 16p '4 al Da 


2 И = 2 
cn(v,, B) HeunG, e E 4h LO 1 3 1, 1 en" (v, f) 
B 165 4 4 2^ 


Рог @(v.)=1 @(K(B))=0=> @(v)- 


On peut également considérer à profit le développement limité autour de la valeur K(6) de 
l'équation différentielle, et obtenir une approximation sous la forme de fonctions de Bessel 1 
lorsque les valeurs propres h sont positives et de valeur importante : 


"x)-4 h+vlv + 2 sn (xB) x)= avec x € lx et > 
у) рна) о ки) а љо 


= Е 5 — => у"(х) | zk dE 
ï =K(B)-x= у") h ! Tue »G)- Ax || ch Bx EK. Liz 


L 
2 


Là encore la solution régulière à l'infini est la fonction de Bessel 1, soit donc l'approximation 
suivante ` y(x)=V%1_,(VhXx)=J/1„(Vhx). Pour la valeur particulière m=0, l'expression approchée est 
Wis 


alors particuliérement simple : 


ss . OE, -0 avee xek, KI е к 
h 


4 cn°( 
T sp (ei). A 
Y= |с» ARTE) 


Construction d'une deuxième catégorie de solution locale réguliére à l'infini par un 


développement de Fróbenius 


La solution locale réguliére à l'infini de l'équation suivante requiert une approche similaire à la 
précédente étant donné que la singularité est également à l'infini : 


PH) + fa @ r) 1 VE avec ue [u,2K(a)] 


Ou sn’ (u,a) 
ü-u«iK(B)  sn(ŭ,a)= : 250 SH Je E kee ml 
a sn(u,a) ou 4 
x-sn'(ia) à 
I = — = sn*(u,a) 
а = = x 


a? 


Ici on peut mettre à profit la transformation homographique a/x qui figure justement parmi celles 
donnant les fonctions de Heun équivalentes, pour obtenir une expression « simplifiée » de la 
solution de Frübenius. Les 2 expressions des solutions 1 et 2 au point singulier infini données par 
R.S.Maier dans « The 192 Solutions of the Heun Equation » en pages 24, 25 et 26 sont les 
suivantes, soit réguliéres soit non réguliéres à l'infini. On remarque que les solutions #1 et 42 sont 
interchangeables par la substitution : a<->B : 


Fonction de Heun point singulier х= oo avec d+ B EDERT EE 


x * HeunG( aq I alala y—ó+1) B + Sba a у+1@-8+1@+80-у-8- :2)- 
Solution #1 — X 


=x нето, aq salto y-ó+1)- Вч DER y *L& — В 4 Lat) 


СО y 8 +1) G8, B.B у+18-@+1,@+8-у ë +1;2)= 


Solution #2 => 


=x ^ HeunGi| a, q+ EMP y —ó4 1) @ + 5} 8,8 y*LlB-à4 Lat) 


Dans le cas d'une équation de Lamé où m=0, les paramètres sont les suivants : 


1 A - o 1 1 
a= = @ Ô =€ 
a? 1 4a? P 4 S 2 
4 en us T aaro 2 


Alors pour m=0, les solutions #1 et #2 sont identiques et toutes deux régulières à l'infini, ce qui 
nous améne à proposer la forme : 


Bain) а) = A> > HE) Le [oe = L kein 


> THW) | fa ü т уни) avec ше[и,.2 (а) 


1 1+a°-4 1 3,1 A 
>H = /sn(u,a)HeunG А —,—,1, sn ,a 
(и) # (2 TE 0 ie ) 


Lorsque m» 0 comme z 1—2" =g et B- 2m +1 o, la seule solution régulière à l'infini est la 
4 4 


solution #2, soit pour l'équation différentielle suivante, le développement de Fróbenius-Heun : 


н) da (m d 1 нидо avec u e[u,,2K(a)] 


ди 4 sn’ (u, 
= 1—2т ~ 2т+1 1 À 1 
a = et B= а = ЙЫ EE 
4 B 4 a? 4 4a? Б 2 


2 2, 
(2т+1) (a LA 2m+1 2m+3 > Zë eina) 


2m+1 
Solution #2 — H(u)= ОДЕ ` 16a° 4 4 


La récurrence du développement autour de l'infini s'écrit alors comme suit lorsque l'on se réduit 
uniquement au cas m=0 : (toujours plus efficace que d'avoir recours aux fonctions de Heun 
implémentées sous Mathematica) : 


CH), ! Á = 
С isto doo eee ика] e HEK()-0 


ди 
н(и)= Joo a )HeunG |2. 1+ a? сш 13,1, E. IB STE Amt bud 
б. Аа о ш Z 
1+a*-42 , RH 
4 =l A= ET 4 Aus 16(л +1) 


2 
немо] r E ihla) 
Pow  H(m)=1 HQK(a))- 0 H(u)- TES oi 16а 


sn(u,,@) 1 1+6?-44 13.1 , 
HeunG j E ‚а 
еип [2 lGa? °g g”? э (u, ) 


Remarque : toutefois, il ne faut pas oublier que si u0 < K(a) alors le développement de Frôbenius 
n'est plus valable car l'intervalle comporte deux points singuliers. Il faut donc dans ce cas trouver 
une solution approchée notamment sous la forme de fonction de Whittaker (voir le paragraphe 


suivant). 


Valeur approchée de la deuxiéme catégorie de solution locale réguliére à l'infini 


On peut également considérer à profit le développement limité autour de la valeur 2K(a) de 
l'équation différentielle, et obtenir une approximation sous la forme de fonctions de Whittaker : 


Tre б r) e Jil avec ue|u2K(a) et H(2K(a))=0 


ди? 4) sn°(u,a 


Notation v >u ц Xo m =À >v(v +1) Hy h<0—1=-h>0 


à H viv +1) x)= avec xelx a 
eich, AED E 


sn*(x,a) (x-2K(a)) = y" (X) fas "Gu, 02 y(x)= Ах Уй Micha JIK | Aï) 


X -2K(a)-x 


Là encore la solution régulière à l'infini est la fonction de Bessel l, soit donc l'approximation 
suivante : (X)- x1 as)- „т, (Vna). Pour la valeur particuliére m=0, l'expression approchée est 
E 


alors particuliërement simple : 


2 


eben == share) < Dou 


P [m 1 -hQK(a)-x)) 
2K(a)-x, AN AQ K(a)- x,)) 
Ce développement n'est valable que si la valeur h en valeur absolue est suffisamment grande, mais 


pas forcément trés importante en soi. Dans ce cas j'ai été surpris de constater que l'approximation 
par les fonctions de Bessel est remarquable sur tout l'intervalle considéré. 


Développement perturbatif pour une approximation des solutions locales à l'infini 


Que ce soit la première ou la deuxième catégorie de solutions locales à l'infini, on a vu que le 
développement au premier ordre met en jeu des fonctions de Bessel l. Cette approximation n'est 
pas adéquate pour des faibles valeurs absolues des valeurs propres h. Je tente donc d'améliorer 
cela par un développement perturbatif. Pour simplifier, il suffit de se placer autour de x=0, ce qui 
correspond à la singularité infini de l'équation de Lamé. Proposons la forme suivante de la 
solution : 


у(к)= (к) (к) fG)g (x) ` o ye) Colle 27698 (к) + f(x)g"(a) 
dall — v(9)- /'(0)g(0)+ /(0)e (0) 0 g(0)= 0 


) 
"posteri aler: 


1 
dont le paramètre de perturbation & serait m° - 4- Soit l'équation différentielle suivante de la 
€ DE, 


АД 
SEN gelt) eee eich Be a- Ale 
onction g(x) : gx; )=1 
една l)e Ur et Veeg(x)-12 g(x)-0 
&,(0)« 0... .. g, [x,)= 0... 


Pour la première équation différentielle de go(x), nous avons : g'(x)- 22 E g'(x)=0= g,(x)- 1, et sans 
méme réfléchir beaucoup, la solution constante est triviale, soit la valeur 1. Pour la seconde 
fonction, l'équation différentielle est un tantinet plus complexe : 


"— du 
go(x)=1 g (x) TH g'(x) Mi [ A 


Comme pour la fonction f(x), développons le second membre autour de la valeur x=0 : 


De méme développons le terme f(x) comportant la fonction de Bessel autour de x=0 : 


ER E hs meet E wl Vu? Til 2m«1 | Qí44,2m*l An, leo 
fe) ine re O и) ug Ws SEU D 


La solution générale de cette équation différentielle est la suivante : „ (x)= 44 e - i e Zei 
m 


Comme la solution ne peut comporter de divergence en x=0, cela implique que B=0, et comme elle 
doit s'annuler au point x=x0, il vient immédiatement : 


EG f> s 
COS een à "fy | 


La solution approchée аи premier огаге de perturbation s'écrit donc : 


i ec . jen avc xef] h>0 e | p. 


ee нео) E 


Pour le problëme complémentaire qui nous occupe, cela donne : 


ZORGC Jo dei? avec xelx,2K(a)] h>0 е Ee 


n (u,a 


_ [2K(a)-x 1,(/h(2K(0)-x)) nn E cs 
x PER Sea | i p Ue) A 


S'il s'agit maintenant de trouver le même type de développement perturbatif au problème suivant : 


6 Le EA A du breit, be ge 


La solution approchée s'écrit alors : 


4) 12 


ui BE M GC Tin RD) 


Fonctions de Lamé-Wangerin 


Les fonctions de Lamé-Wangerin ont été introduites pour résoudre des problèmes aux limites dans 
systèmes de coordonnées orthogonales basées sur des surfaces cyclides, telles qu'elles avaient été 
introduites en tout premier lieu par Albert Wangerin en 1875. Une étude de la séparation de 
l'équation dans un tel système de coordonnées cyclidiques, comme celui d'E.G.C.Poole conduit à 
deux équations de Lamé séparées dont l'une a comme particularité que la variable peut être mise 
sous forme purement imaginaire. De ce fait la solution du problème est susceptible d'atteindre la 
singularité à l'infini de l'équation de Lamé sous forme algébrique. 


Dans ce cas une étude des contraintes de continuité de la solution conduisent à établir la nécessité 
pour une des solutions de l'équation de Lamé de vérifier la propriété suivante : 


Remarque k>a k'=\1-k f 


DO ey) fini Lorsque z=+iK(B) 


Pour cela plusieurs études sont dues à Erdelyi dans l'ouvrage « HIGHER TRANSCENDENTAL 
FUNCTIONS VOL Ill > dans le chapitre XV sur les fonctions de Lamé, dans lequel il propose un 
développement des solutions de Lamé-Wangerin lorsque l'argument est principalement de la 
forme z=x+iK(6). Les développements sont alors donné sous la forme : 


Remarque k—a К В 


e (eren Ea 22] 


en(z,a)-isn(z,a) "Cen(z,a)-isn(z,a) 


г=0 


ЯАР e ЕЕ Sr (е6 (е а Ё T eJ sr (tease) 


n(z, a)- isn(z, a) cn(z,a)-isn(z,a) cn(z,a)-isn(z,a) z0 cn(z,a)-isn(z,a) 


Par les équivalences d'expressions suivantes et l'introduction de la fonction amplitude de Jacobi, on 
peut déjà simplifier la forme des développement : 


LR. en(z,a)= Sin(£) en(z,a)+isn(z,a) _ Sin(£)+iCos(£) _ _ Cos(£)-iSin(£) _ äu 
= ES Ge: E = Cos(¿ ) en(z,a)-isn(z,a) Sin(£)-iCos(£) ` Cos(£)+ iSin(Z ) ° 
cn(z,œ)-i B sn(z,k) en(z,a)+i B sn(z,a) _ Jen? (z,a)+ B sn'(z,a) "се +isn(z,a) z 
TERES | en(z,a)-isn(z, a) cn(z, a)-isn(z, a) en(z,a)-isn(z,a) yan п(2,а) —isn(z (2,а) 


FE?" (z, a) zi SA (- ly pem ег (v+2r+1) 2 A, е“ (v+2r+1) 
r=0 
Eae dn(z,a) Y B (- i^o uA aon = dn(z,œ) Beier 


r=0 r=0 


Voyons si ce développement se simplifie lorsque l'on applique le changement de variable z=ix. Il 


s'avère être réel (à une constante près) lorsque l'argument de la fonction est purement imaginaire. 
A savoir que l'on peut écrire : 


sn(ix,a)=i A cn(i x, k)- I dn(i x,k) = dn(x, B) Ballo" 


cn(x, B cn(x, B) cn(x, B) 
1 
en(ix,@) _ _ i Ao аа) а), AA 
sn(ix,a)  sn(x,B) enfix,a)-isn(ix,æ) l a Í+ sn(x, B) 
sn(x, B 
2 К 7 _ dn(x, B)- о sn(x, B) _ cn(x, B) 
dn(i x,a)+i a sn(i x, z) HP) Ab) ET) 


Et qu'à partir des deux développements des fonctions de Lamé-Wangerin, on arrive aux fonctions 
suivantes : 


F?” Guje S geite F2"*(z,a)= dn(z,a) Y` B e rx ce 5 am(z,a) 
r=0 r=0 


(v+2r+1) 


F” (ix,a)= DE LES 


z-ix Ç =7-amix,œ) e = КЕ en 8 * mu (v+2r+2) 
OUS Fx Ë CB) y _ 1—sn(x, B) 
М ї en(x,B) 2 "V 1+sn(x, B) 
(v +1) 
T z i C22 a _ (1=sn(x, B)) ° < „(1 п(х, B)Y 
F"Gsa)-e F, EES 0 41) > 
= (v+2) dn(x, B) ( 1— sn TES 1—зл(х, B) Y 
T2m+ = > 2т+1(. "Ax, TAX: жч —1) PAR 
F (x,a)=e F; (i x,a)- SD sn(x, ed > В,( 1) F 


Et l'on peut affirmer que les deux fonctions sont solutions de l'équation différentielle de Lamé 
suivantes : 


v+2r+l 


2007022 eL) а-о» жа deen | 


en (x, B) 


v+2r+2 


S2m+1 _ dn(x, p) , 1—- sn(x, B) 2 
EU e) SE 2,5 y| Sen 


" ; . 2m (x,a) pru a) 
Condition de Lamé Wangerin -* ——- et — 


Jen(x, В) dent). continue en x — +K(B) 


Toutefois pour l'instant je doute de l'utilité pratique de cette construction dans la résolution de 
problèmes aux limites pour le changement de variable z=ix, sachant que certes les fonctions 
possédes les bonnes propriétés en x=K(6), mais elles divergent en x-K(6). D'autre part je ne peux 
construire de fonctions paires ou ur sans rencontrer la méme divergence. Les valeurs en x=0 


sont les suivantes ` F?"(0,@) = YA F= (0, a) - XB . Et les dérivées premières ne 


s'annulent pas plus en x=0. En x-K(6) nous avons : 


Li, € (s aere air Era» Ean E] б-ки)" = К”(К(д}о)-о 
= Era) Zar äi (к-к Enea 


r=0 


0 


Il est donc clair que cette construction ne peut représenter de solution d'un probléme de Sturm- 
Liouville car les fonctions n'ont pas de conditions aux limites homogène sur l'intervalle [0,K(8)]. 


Et si nous appliquons un changement de variable de la forme z=x+iK(8) comme initialement dans 
les paragraphes consacrés aux fonctions de Lamé-Wangerin par Erdelyi dans l'ouvrage « HIGHER 
TRANSCENDENTAL FUNCTIONS VOL III » , chapitre XV sur les fonctions de Lamé, alors nous arrivons 
aux expressions et aux fonctions respectant l'équation : 


(Кво). exxeik(Ba)s 269). cur a) - 20099) KE zk Е 


a ѕп(х,а) ia ѕп(х,о) іѕп(х,а) n'(x,a) 
Avec ce changement de variable z=x+iK(6), les fonctions construites ont maintenant les expressions 
suivantes : 


a sn(x,œ) ` 1— dn(x, a 
1-- ап(х, a) ENS e 


do. v+2r+l 
rel E да v+2r+l sn(x, a) A, - dn(x, a N + 
(s a) 24 | ) -È 1+ pe ) 


e = -i(en(x- iK(B) o) i sn(x+iK(B),œ)) = 


1+ ап(х,а 
E v-2r«l MEER EXE 
р?! T Р a)Y Ba"? (sn(x, о)) E en(x,a) + (aea) МЕ «сп(х,о) ўв сања) 2 
I і Ки s (1+ ап(х,а)) "` — sn(x, a) £x "Ae dn(x,a) 1+dn(x,a) 2 \1+dn(x,a) 


E (x, a) di p (x, a) 


Jsn(x,a) Jsn(x, (x, a) 


Et correspondent aux deux problémes de Sturm-Liouville suivants : 


DE e) CE E et y (K(a))-0— F2"(x,a) 


sn*(x,œ) у(0)=0 er у(К(а)) =0 E?" xa) 


Condition de Lamé Wangerin continue en x = 0 


La récurrence des coefficients du développement de ces fonctions est la suivante (formules 15 et 
16, chapitre XV, < HIGHER TRANSCENDENTAL FUNCTIONS VOL III >) : 


Н -2h-v(v +1)a? 

4,20 

"i (+1) (2-0? ))4, «Qv +3)a24, =0 

(2r -1)v +r)a? EE (2-a?))4, +(r+1X2v e 2r-3)o? A 
B,-0 

(H - (v +2} (2-a2))8, + (2v - 3)? B, = 0 
(2r+1)v+r+)a?B, (н - (v 2r 25 - o?)B, +(r +1)2v +2r+3)a28, = 0 


=0 


rl 


Constructions de H.Volkmer dans < Eigenvalue Problems for Lamé's Differential Equation », 2018 


H.Volkmer note dans son article : 


lëfteg ue geng ale зы sns 
en(z,a)-isn(z,a) 

1- 1+ 4 
ъ= ЮЛ т = 1— bed nm =1 ВЕ EE 


a? i^ a^ lf, 
Р 0, > + 1 +2n +1)-1 
our n e [0,m.] >= É + alt п+1) 


а? 1 а? 1 4 а? 1 
De plus (m -nXn. -n)= lem a — +27 |= (oem +2n)-1 
4 n 4 n a 4 n 


De plus nous avons : 


Jn = sn(z,a)-icn(z,@) —= = sn(z,a)+ ісп(2,а) 


2 
= үп + Е = sn(z, a) => a?sn'(z,a)-12 -ат?(2,о)= € d n), =n) 
Jn n 
> idn(z,a) "тут nn 


La forme trigonométrique de l'équation de Lamé à l'aide de l'amplitude de Jacobi s'écrit : 


CE 2 = am(z,œ) = Cos(£ ) = sn(z,k) 
(i-a? cos (С) (c) o Cos(£ )sin(£ )v (£ )+ (n -v(v + )о°Соѕ?(6)) (С) = 0 
e irao [eco OL GO) tlc" )=0 


Avec la variable n, l'équation de Lamé devient pour diverses formes de la fonction y(n) : 


n s Mn JE. с E Ae ge 


dy — 2m n-m n-m) dn 41 
vo) cien n ЕШ A L 4 A AND party E, te 
y(n)= Jn, -ng(n)=> onm =m. Yn rh Aem. cvy jeta E (n dai? 


(n)= Jm nm =ng(n)=30r a(n +n,)=22-a?) 


o^n(q-n Xn dE JA DÉI var i m T) tn 


dy  2im m-m n-m) dn 


Een ЗЕЕ ВЕ ER air, eh Ca 


Avec le changement de variable z=x+iK(6): 


Г. р Ud _ dn(x,œ) 
sn(x K(B), ) ia sn(x,a) 


а ѕп(х,а) 


сп(х+іК(В),оа) = 


Nous pouvons ёсгіге les deux expressions suivantes : 


CRU a eee ky) e * este aenea) 


ant + iK(B).a)- en(x + iK(B), œ) = SE еее сп(х,а) 


H.Volkmer donne раг ailleurs les deux identités suivantes : 


q-aqm-n4m-"n Үп 0<n<m 
VJ1-n +q 41-75-74 = 24A - Bm. =n 
J1-n+q - 4172-4 = A2 V B n, =n 


Pour s'en convaincre il suffit de multiplier les deux identités entre-elles, et de constater que la 
deuxième identités s'annule en n=n1, la factorisation devient alors évidente. Comme : 


E а NTh Ny To” . 1 1 
idn(z,œ)=> n Sg CS et EE) 2 


1 venil те enne) 


1| 1 А 
1-7 all L DE 2i dn (dn(z, o) en(z, a) 
m E ace 


iz 1 
1-7 -q= Узе *n* Jdn(z,a)- cn(z,a) 


Nous arrivons aux identités suivantes : 


i E Jdn(z, a.) en(z, a.) Ej Jdn(z,a)) - en(z, a) 
i MES 


M n = 


Jm =n x j* ei Van(z,a)+ ME Jdn(z,a) - en(z,a) 


Notons une petite erreur de calcul de H.Volkmer dans son article de 2018, « Eigenvalue Problems 
for Lamé's Differential Equation », en page 11 pour les deux formules 4.7 et 4.8 de signe centrale à 
inverser. 


Soit l'évaluation pour la valeur de l'argument z=x+iK(6) pour lesquelles les expressions sont bien 
réelles comme attendues : 


Si 2=х+1К(В)= Jm =n р! Т 
i Van(x,a)+ a сп(х,о) + J dn(x,a)-a en(x, a) | 1-dn(x,a) оп? (х,а) 


Jh B Ja sn(x, a) die 1+dn(x,a) (1 dn(x, a) 


= Ji —-T =n 


Idem pour {= 


| Jdn(x,a)+ a сп(х,о) — J dn(x,a)- a cn(x,œ) 
_ YEA Лев. aqhapa) 
=== 2 J dn(x,a)+ a en(x,a)+/dn(x,a)- a cn(x,œ) 
Er J1- B J1+ ап(х,а) 


Les deux nouvelles constructions que propose H.Volkmer sont les suivantes : 


r=+œ viir) 


e. F” (z,a)- n -n YX. An ° 
B dE ansa) r-0 


r=+0 v+2r+1 


Е2"1(2,а)= {п-п Y. Bn ? 
r=0 


v+2r+1 


Ets a) = J dn(x,a)+ a en(x,o) + Jdn(x, a.) - a сп(х,о) 5 en 2 
. Geck J1- B [I+ an(x,a) per 
z=x+iK(B)= ENT 
or _ Jdn(x, a) a en(x,œ) =J dn(x,œ)-œ en(x, œ) 3 i-es 2 
ge (ъа) + B J1+an(x, a) 25 
Persa) Fa) 


ti =0 
BE e Ба) continue еп х 


Condition de Гатё Wangerin 


Avec pour récurrence des coefficients les expressions suivantes : 
H=2h-v(v+l)a? | B-41-a? 
1 2 
A,=0 е E (2 WE | (v | : | | 2p(v+2) A,+(2v +3)а?А =0 


2 
etra [н (2 eJ ( 2r4 d | -25 ЕЕЕ =0 


4 
2 
B,=0 а [ (2 ejes | asv] B, - (2v +3)a?B, 20 


2 
оез дазв (н (2 eJ ( 2r4 d | 2B| v +2r+ 312 +(r+1X2v+2r+3)a?B,, 20 


4 


Il est à noter que la matrice de récurrence est symétrique. 


Les premiëres constructions dont parle H.Volkmer sont équivalentes qux constructions proposées 
par Erdelyi dans l'ouvrage < HIGHER TRANSCENDENTAL FUNCTIONS VOL III > , chapitre XV sur les 
fonctions de Lamé 


r=+0 у+2ғ+1 


F?"(2,a)= 241 ° 


п=е 
p Mox "ув n = 


v+2r+l 


еа) Sra (cien) 7 
Z= x+iK(B )= Vn 7 TD п = 2065 2 2) E B Kë у+2г+1 
+ dn (x ol _ e 
POmH cn(x,œ) ZE 1-dn(x,œ) 2 
F; (x,a) = B, 
1+dn(x,œ) Pari 1+dn(x,a) 


posa (x, a) 


sn(x,a Jsn(x,a ) 


Condition de Lamé Wangerin continue en x = 0 


Je rappelle la récurrence des coefficients du développement de ces fonctions est la suivante 
(formules 15 et 16, chapitre XV, « HIGHER TRANSCENDENTAL FUNCTIONS VOL III ») : 


H=2h-v(v+l)a? В=ҹ1-а? 

A,=0 

(н - (v +1) el, - Qv -3)a? 4, =0 

Qr v e r)a?4, +H (и +1+2r)(2-a2))A, 12v 92r +3)a24,, =0 
B,-0 

(H - (v + 22 -a?)B, + (2v - 3)? B, = 0 

Qr v r1), ,+(H -(v 2r 25 -o?)B, +(r+1)(2v +2r+3)a2B,„ = 0 


Passons à la première construction des fonctions de Lamé-Wangerin toujours proposée par Erdelyi 
dans l'ouvrage < HIGHER TRANSCENDENTAL FUNCTIONS VOL III > au chapitre XV, sous la forme : 


r=+œ ege 


p (z,a) = Y. 4A, (sn(z,a.)) (2) = cn(z,a) Y` C, (sn(z, oy (202) 
r=0 0 


F2" (z a)= dn(z,a) Y` B, (sn(z, a) A _ dn(z, a en(z, a) Y` D, (sn(z,œ)) (27) 


F2"(2,0)= eR S4 m. a)" = RU Ca)” 


r=+00 


Elte ll Sema)" анана | ] Lora)” 


Pour un argument purement imaginaire z=ix, il vient : 


j _ sn(x, B) nli x, k) = n(i x _ а(х, B) 
nuu) cn( ,k) EXE dn( ,k) Aa; В) 


) 
O m Z бй = AR " 
FP(ixa)e ? nr Ўл) |218) _ xo спі В)" So (1 (een 


ta 
S 
~ 
x 
`œ 
x 


E d E EE 


sn(x, ) sn'? (x, B) par; 


sn"? (x, В) < r 


Р'"(у ке 2 Fix Е en Qe By, Ë у cn(x, ) 2r 
i í i Á í ü sn" (x, B) r=0 ' 


Sp iz (v+2) an T ” 
Раа) E" (xa) MES ED F p y| ten) 


Même en obtenant une construction, de valeur réelle, qui respecte la condition imposée cette fois 
en x=K(6) et x=-K(6), il y a divergence en x=0. 


Pour un argument de la forme z=x+iK(6), nous obtenons : 


marria- enfrrik(pha)= Së). ap) Í - Sal 
a sn(x,a) іа sn(x,a) isn(x,a) 
E?" (x, a) - E?" (x iK(B).a)- S Aa? sin sar) 
r=0 
F2? (x, a) - iE?" (x + iK(B) a) = en(x, a) Y` Ba" sn"? (x, a) 
r=0 
ou encore selon la deuxiëme forme proposée par Erdelyi : 
Fer (x, a) = iF?" (x +iR(Bha)= dal el Y. C, a?" sn?" (x, a) 
r=0 


F2" (x,a)= E?" (x +iK(B) a) - –сп(х,а)ап(х, a) Y D arteria rare (х,а) 


r=0 


Et correspondent comme avec la forme de développement trigonométriques aux deux problèmes 
de Sturm-Liouville suivante : 


sn (x, a »(0)=0 et y(K(e)=0— А" (x,a) 


DE ED e= RT et y'(K(a))2 0  F?"(x,a) 


Voici les relations de récurrence des coefficients Ar, Br, Cr et Dr : 


A,=0 
(n-(v +1) (1+a2))4, - 20v +3)a°4 =0 

(v+2r-1)v+2r)4, (а - (v e 271 (а?) 267v + 2r + 3)a°4 
B,-0 

(a- (v +2) -(v frech, +2(2v - 3? B, =0 

(v+2r)v+2r+1)B, (8 - (v 272 - (v «2r 1 a2)B, - Qr 12v +2r 3) B. 
C,-0 

(a- (v 1) -(v +2Y a Ju + 2(2v +3)aC, = =0 

(v +2r)(v +2r -1)C,. OCH (v+2r+1) y -(v+2r+2) a*)C, «(2r Av +2r +3)a?C,„ = 
D,=0 

(a- (v +2) (L2), + 2(2v +3)a2D, =0 


(v 2r 1v -2r2)D, +(й—(у 2r 2) (i +2(r+1)2v +2r +3)a2D 


+1 


=0 


r+l 


=0 


De ce fait on se trouve donc un peu démuni avec le changement de variable z=ix (problème aux 
limites dans le système « Flat-Ring » d'E.G.C Poole) , car aucune méthode de construction valide ne 
semble émerger de la littérature. 


Sauf peut-être dans l'article du même Erdélyi de 1948, < Lamé-Wangerin Functions > dans lequel ce 
dernier propose cette fois pour le problëme électrostatique dans le systëme de coordonnées < Flat 
Ring > d'EGC.Poole, les développements suivants : 


е ( dn(z,a)-aœ en(z d 
Е2"(2,а) | и 
E B 


os AY dn(z,a)- a сп(2,а) У 
Е; (z,a) 2 sn(z,a) В, B 


Ce qui donne pour un argument purement imaginaire : 


.sn(x, ) | a 1 | = dn(x, p) 
sn n(i x, a)=i at ) cn(i x,k) = SESCH dn(i x,k)= cn(x, В) 
2т(. zm, E À, dn(x, B)- а E bed TE v+2r+1 cn(x, gr 
F, (i x,a) > jum | cn(x, В) J 2,45 (dn(x, flot" 


ysadi v+2r+2 PEN v42r4l 
pu А j E sn(x, B) B, (ze "d E | В. у+2г+2 cn(x, p) 
v (ix a) " en(x, В) Z В??? cn(x, B) i sn(x WK „B (dn(x, B)+ a) ^"? 
F2" (x a)= EP (ix, a) )= 25 В?! en(x, B)?" 
v RN (dn(x, B)+ a) ^ 


poo v+2r+l 
pou : = rg" ok : B В" cn(x, 8) 
NE CD 


La construction semble présenter les bonnes propriétés de finitude sur l'intervalle [-K(6),K(8)] et à 
ses bornes : 


m sn(x, В) ы p2m z NS v+2r+1 cn(x, py: oc ve 
у> 7 et EES La (x, a)= BI x NAE 24 p" (ано В) а)" сп(х, В) 

_1 sn(x, p) : T'2m+l 2 AG v+2r+2 сп(х, ву oc ve 
v2 2 et cn(x, B) F; (x, a)- |sn(x, B): sn(x, B) =з В,В (dn(x, B)- a) ^? сп(х, В) 


Pour v>-1/2, alors les deux fonctions s'annulent en x--K(8) et x-4K(8), et pour v=-1/2, la valeur 
reste fini (mais cela ne garantie pas qu'il y ait continuité lorsque l'on passe de x=-K(6) à x=+K(8)). 
Les valeurs en x-0 sont également de nature homogénes puisque les développements sont paires 
ou impaires: jg2"(9,5)-9 ` g?""(9,4)-0- Soit en récapitulant pour les deux bornes de 


l'intervalle : 


F?" (0, a)= 0 Ке (0,а)= 0 
Е2"(К(8),а)= 0 FE" (K(B)a)-0 
Ё?" (=x х,а a)= Ё?" (х,а) Tn (= x,a)= -F?""(x.a) 


Il y a donc de bonne chance que cette construction puisse être utilisée dans un problème aux 
limites de Sturm-Liouville, comme celui sur de l'iso-surface u=Cste portée à un potentiel 
électrostatique constant. 


La récurrence à appliquer pour la détermination des coefficients des deux développements : 


Ein (x a)= Уа Bree сп(х, py 
v > Z r (dn(x, Glaf" 


F2m+1 у+2г+2 сп(х, Bl 
Е, (x, a)= s n(x, B DE „В (ап(х, B)-- a) ^ 


(vv +1) (и +1) )1+a2)-2h)4, +(2v +3)B24 = 0 
Qr 1v e 824, ,+(v(v+1)-(v +1+2r) fi ect 2)4, (o 1v +21 +3)824,„ =0 


B,=0 
(vv -)- (v 2 fi a2)- 22)B, - Qv +3)828, = 0 
Qr 1v +r+1)B*B, о Y1+a2)-2h)B, +(r+1)2v+2r+3)828,., =0 


Pour obtenir cette récurrence, il suffit d'injecter les deux formes : 


rre) S a (sesion 
r=0 


r=+œ v+2r+2 
E a) sn(z,a) N „ [are еа!) 
r=0 


Dans l'équation de Lamé originale, et d'identifier les puissances de la variable „ = dn(z,a)-a cn(z,a) 


terme à terme. Pour cela il est particuliérement utile de rappeler les correspondances de variables 
suivantes : 


Ç —-sn'(za) et u=dn(z,a)-a cn(z,a)- J1- o? - 4l et В= 41-а? 
Soit Ç =1 s É u? f =>. -2 et er 
u 


en(z,a)= Al £ e 
ап(2,а)= J1-a?£ te 


=? 2 2 
sn(z,a)- Jc = 411 (6 e) 
4u'a 


E eter ica E 
u 


Là encore on peut affirmer que ces deux solutions respectent bien l'équation de Lamé modifiée : 


Revenons sur la propriété de finitude imposée aux fonctions de Lamé-Wangerin, sous la forme : 


1 
Isn(z, a? y(z) fini Lorsque z =>iK(B) S Lim ЕЕ: ОИ 


хә+К(8) [cnlx, B ) 


que nous venons de vérifier. Une fois injecté la forme de la fonction dans la solution complëte de 
l'équation de Laplace, il faut tenir compte du terme de séparation. Voyons juste le comportement 
de la solution en v autour des valeurs v=+K(6) et v=-K(6) : 


l убу) 
Т R —>T —— 
x2 ( ) (vv) ) Гоу, B) =? (v)œ e. P) 
у(у)= E” (v, a)- S'y p cn(v, B^ i au cn(v „BY ñ: 
LI r (dn (v, B )+ Sg ` 0 (2a y +1 
xtK(B) ou 
r=+o% у+2г+1 vH 
E pro А В v+2r+2 en(v, B) = +B v+2 спі, B) 
LOE (v, a)= sn(v p) m „В (dn(v, В)+о) „8 Ray” 
1 
0 si v>-- S 0 si v»-— 
2m 2т+1 
> Li kd =] 4 gr y à Lim Eva) =) p gre Q | 
хэ+К(В) [ед v, B „8 “ы x5K(5) dent. B $ „B j v2-— 
Qa)" 2 (a) 2 


La solution complëte avec la fonction impaire présente un léger inconvénient de discontinuité entre 
les valeurs v=+K(6) et v=-K(6). Toutefois cette discontinuité est compensée par l'annulation de la 
fonction de la deuxiëme variable sur le plan équatorial, comme nous le voyons lors de la résolution 
du problëme aux limites électrostatique (voir ci-avant dans la section < Systëme de coordonnées 
d'E.G.C Poole). 


